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SUMMARY 


This report discusses activities under NASA Contract NGR 39-009- 
077, involving a theoretical and experimental study of high pressure 
fuel droplet combustion. The investigation was divided into two phases, 
as follows: 

1. Droplet combustion at elevated pressures (’-100 atm) in air 
maintained at room temperature. The droplet was in a natural 
convection environment with Grushoff numbers in the range 10 4 - 
10® (based on approach conditions.) 

2. Droplet evaporation and combustion at elevated pressures (1-40 
atm) in a high temperature (600-1500° K) flowing combustion 
gas environment which simulates actual combustion chamber con- 
ditions. The ambient oxygen concentration of the droplet was 
varied in the range 0-13% (molar) with forced convection eval- 
uated over the Reynolds number range 10-800 (based on approach 
conditions) . 

Porous spheres were used to simulate fuel droplets in order to de- 
termine burning rates and liquid surface temperatures under steady con- 
ditions. The fuels considered in the tests included methanol, ethanol, 
propanol-1, n-pentane, n-heptane and n-decane. 

The theory used for comparison with the experimental results was 
based upon an extension of the variable property, variable specific heat 
model of droplet combustion proposed by Goldsmith and Penner. This ap- 
proach has the advantage of postulating realistic property variations 
for the fuels considered in the tests. Constant specific heat and con- 
stant property versions of the theory were also compared with the ex- 
perimental results. The effect of natural or forced convection was 
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treated by means of a multiplicative correction based upon the Gras- 
hoff or Reynolds number evaluated at approach conditions. Phase 
equilibrium at the liquid surface was determined by both the conven- 
tional approach typically used at low pressures, as well as a high 
pressure version which allows for phenomena which become important as 
the thermodynamic critical point of the liquid phase is approached 
(solubility of ambient gases in the liquid phase, etc.). The phase 
equilibrium calculations of the high pressure theory employed a modi- 
fied Redlick-Kwong equation of state with mixing rules for multicom- 
ponent mixtures. 

The conclusions of the study are as follows: 

1. Both the low and high pressure versions of the theories gave 
essentially the same burning rate predictions over the range 
of pressures for which they overlap. The greatest difference 
between the theories is that the low pressure model predicts 

a significantly lower pressure foi he liquid surface to reach 
its thermodynamic critical point for given ambient conditions. 

2. The various property treatments; variable property-variable 
specific heat, variable property-constant specific heat and 
constant property models gave similar predictions as long as 
any assumed constant property was evaluated at an average con- 
dition. 

2. Discrepancies between theoretical and experimental burning 

rates were similar to those encountered for atmospheric pressure 
tests. The theories modelled the effect of variations in tem- 
perature, pressure, ambient cxygen concentration, droplet dia- 

pri-eeuunG page blank not filmed 
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meter and convection conditions reasonably well. The repre- 
sentation of the effect of fuel type was less satisfactory 
with the theory giving good predictions in some cases and 
errors up to 50% in other instances. 

4. At pressures up to 20 atm, liquid surface temperatures were in 
good agreement with the low pressure version of the theory. 

At higher pressures the surface temperatures tended to approach 
the surface temperature predictions of the high pressure theory 
in most instances. 

5. The experiments gave evidence of the liquid surface reaching 
its thermodynamic critical point for the combustion of methanol 
and ethanol at pressures in the range 80-100 atm. This pres- 
sure range was in agreement with the critical combustion pre- 
dictions of both theories. 

6. A major advantage of the high pressure model is that it pro- 
vides a more satisfying physical description of conditions at 
the liquid surface as the thermodynamic critical point is 
approached. With the high pressure model, the composition of 
the gas and liquid phases approach one another with increasing 
pressure, becoming identical at the critical point. This pro- 
vides a smooth transition into the compressed gas combustion 
regime at high ambient pressures, rather than the artificial 
composition discontinuity that is implicit in the low pressure 
model near the critical combustion condition. 



CHAPTER I 


INTRODUCTION 

1. 1 General Statement of the Problem 

Many important combustion devices such as the combustion 
chamber of diesel engines, gas turbines, liquid propellant rocket 
engines, and industrial furnaces involve the introduction of the liquid 
fuel into the combustion space in the form of a spray. The 
heterogeneous combustion of a liquid fuel spray is a complex phenomenon 
and a realistic description of spray combustion requires a knowledge 
of the mechanism of evaporation and combustion of individual fuel 
droplets. Therefore, an understanding of the behavior of fuel 
droplets in a hot, oxidizing gas is useful to the rational design of 
efficient combustion systems for these applications. 

The combustion of droplets has been the subject of 
numerous investigations, with particular emphasis in recent years on 
high-pressure combustion. High pressure combustion, especially 
in the vicinity of the critical point of the fuel, gives rise to a 
number of problems not encountered with processes at low pressures. If 
the critical temperature of the droplet is approached, phenomena such 
as fuel stripping from the drop may become important due to reduced 
surface tension. If stripping is absent, the droplet becomes a puff of 
gas when the critical condition is exceeded and burns in an unsteady 
diffusion flame. Any of these instances give rise to conditions where 
the droplet combustion process is no longer vaporization controlled, 


as at low pressures. 
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A number of models have been developed to describe the 
combustion process at high pressures, but the experimental 
investigations to date have not considered a wide enough range of 
conditions to provide a stringent test of high pressure theoretical 
predictions . 

In addition, most of the high pressure research to date has 
focussed on a quasi-steady analysis of droplet combustion in 
essentially a quiescent gas environment. In studies that considered 
the effect of convection, the ambient gas temperatures were generally 
not high enough to realistically simulate actual high pressure 
combustion. More importantly, previous studies of high pressure 
droplet combustion did not provide steady state burning rate data, 
giving rise to uncertainties in the interpretation of the data. The 
present investigation of high pressure droplet combustion under free 
and forced convection conditions was undertaken to provide fundamental 
data and further insight into this problem, which is of great 
practical importance. The investigation concentrated on temperature 
and pressure levels high enough to simulate an actual high pressure 
combustion process. 

1. 2 Previous Related Studies 

1.2.1 Experimental Techniques in Droplet Vaporization 

Past experimental investigations on the vaporization and 
combustion of liquids have employed three basic techniques. These 
include the suspended droplet technique, falling droplet technique, 
and porous sphere technique. 
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The suspended droplet technique may be used to obtain the 
rate of change of droplet diameter or size as a function of time. The 
fuel droplets are usually suspended on stationary quartz fibers and 
the rate of combustion after ignition is observed photographically. 

The suspended droplet technique yields accurate results for the 
variation of the droplet diameter as a function of time. Some of the 
earliest experimental investigations of droplet combustion utilized 
this technique. [1, 2, 3, 4] 

In the falling drop method, small freely falling fuel 
droplets are ignited and observed during burning, usually through 
photographic methods. This technique is particularly suitable for 
studies of very small droplets. In practice a single droplet or low 
density cloud of droplets is produced by a suitable generator, such as 
an electrostatic generator [5], ultrasonic atomiser, vibrating steel 
tube [6], or spinning disc atomiser [2, 7]. Usually, the droplets are 
allowed to fall under gravitational forces or are projected into a 
suitable hot environment where self-ignition occurs. 

The porous or supporting sphere technique provides a method of 
studying the steady-state combustion of simulated droplet burning. The 
supporting sphere technique consists of supplying liquid fuel to the 
surface of a supporting, inert, solid sphere at a rate equal to its 
iT?te of combustion. Spalding [8] was one of the first investigators 
to use this method. In porous sphere studies the fuel is supplied 
internally to a stationary, non-reactive, porous sphere at a steady 
rate which is just sufficient to maintain a liquid film on the surface 
of the sphere during burning. 
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The porous sphere method is convenient in that different 
diameter support spheres may be used for a variety of experiments 
involving the steady-state combustion of fuel droplets. Most 
importantly, of the three basic droplet burning techniques only the 
porous sphere method is a truly steady-state process. 

1.2.2 Droplet Evaporation 
1, 2.2,1 Low Pressure Investigations 

Previous investigations on pure droplet evaporation mainly 
employed the free droplet and porous sphere techniques. Ingebo [9] 
considered nine pure liquids and a range of air-stream pressure, 
temperature, and velocity conditions in evaporation studies. He found 
that the droplet surface temperatures were approximately equal to the 
wet-bulb temperatures . 

Investigations by Ranz and Marshall [10] confirmed the analogy 
between heat and mass transfer at low Reynolds number, and verified the 
simple expression for the Nusselt number at zero Reynolds number. Based 
upon their experimental findings, Ranz and Marshall developed a general 
correlation for the vaporization of spherical particles. 

Studies by Borman, at al., [11] considered hexane, decane, and 
hexadecane drops vaporizing in air at various temperatures . The 
pressure, range was from one to five atmospheres. Their experimental 
measurements of droplet temperatures and mass transfer rates compared 
favorably with theoretical predictions over the pressure and temperature 


range investigated. 
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1. 2.2.2 High Pressure Investigations 

Torda and Matlosz [12] performed a series of high-pressure 
vaporization experiments for n-pentane droplets. In these experiments 
droplets were injected into and supported in a heated and pressurized 
test chamber containing nitrogen. The pressure range for the data was 
200-1400 psia at a fixed temperature of 390°F. The flow field 
environment was essentially due to free convection, The attempts of 
these investigators to compare their theory with the experimental 
results were not particularly successful. These experiments involved 
pressures that were in excess of the critical pressures of the pure 
fuels under test. 

The investigations by Savery and Borman [13] considered the 
histories of vaporizing n-heptane and Freon-13 droplets suspended in a 
heated air stream. The temperature range for their tests was between 
100°F and 300°F and the pressure range was from 1.5 to 100 atmospheres 
Comparisons were made between the measured equilibrium temperatures 
and vaporization rates, and analogous values predicted from a quasi- 
steady theory uncorrected for high pressure properties. A corrected 
version which included the effects of ambient gas solubility and 
other high pressure corrections to the thermodynamic properties was 
also considered in the comparison. Some of the conclusions of this 
investigation are: 

1. The corrected theory gave good agreement with 

measured equilibrium temperatures and predicted 
vaporization rates (within 25% of the measured values) 
at low ambient temperatures. 



6 


2. The uncorrected theory gave good predictions of droplet 
equilibrium temperatures but underestimated the 
vaporization rates by as much as 80% at high pressures. 

3. The mass transfer rates predicted by the corrected 
theory agreed tc within 33% of the measured value at 
reduced pressures below 1.5. The predicted mass 
transfer rates become progressively lower as the 
pressure increases. 

4. At moderate ambient gas temperatures and high pressures, 
the corrected and uncorrected theory predictions were 
lower by 1G-15°F than the measured equilibrium 
temperatures . 

More recently, Tarifa [14] considered theoretically the 
deviations from the quasi-steady theories for subcritical and 
supercritical droplet vaporization and combustion. He concludes that 
the corrections to the quasi-steady theory are of the order of the 
square root of the gas density to the liquid density. 

Matlosz, Leipziger, and Torda [15] conducted an experimental 
and analytical study of the evaporation of n-hexane droplets in a 
nitrogen and argon gas environments. The ambient gas temperature was 
540°K and the pressure range was from 6.8 to 102 atmospheres. In 
agreement with other investigators they found that the effects of non- 
ideal behavior of the gas phase was important at high pressures. 

In a theoretical analysis of the vaporization (with and 
without combustion) of dodecane droplets in heated air at high 
pressures, Rosner [16] demonstrated certain inadequacies of the quasi- 
steady approximation. He concluded that the evaporation process for 
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isolated fuel droplet just below its critical temperature cannot be 
reasonably treated as quasi-steady even in cases where the boundary 
conditions for temperatures and compositions are strictly steady. 

1.2.3 Combustion Under Natural Convection Conditions 
1.2. 3.1 Low Pressure Investigation s 

The earliest droplet combustion studies concentrated on the 
steady burning period of the droplet. Analytical studies by Godsave, 
Spalding, Hottel, et al. , [1, 7, 8] gave predictions of burning rates 
that were in reasonable agreement with experimental results conducted 
at atmospheric pressure. [2, 3, 4, 7, 8] These studies were based upon 
a quasi-steady analysis that assumed: 

1. Steady state ccnditions for fixed droplet sizes. 

2. Temperature and concentration profiles in the gas phase 
adjust instantaneously to changes in the boundary 
condition at the droplet surface. 

3. The regression rate or radial interface velocity of the 
droplet surface is small compared to che velocity of the 
vapor leaving the liquid surface. 

4. Reaction rates at the flame zone are fast compared to 
the rates at which reactants are transported to the 
reaction zone. 

5. Fuel and oxidizer combine in stoichiometric proportions 
at the flame surface. 

In Godsave* s analysis the temperature dependence of the 
physical properties was neglected and constant average values were 
employed. Goldsmith and Penner [17] eliminated some of the restrictive 
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assump lions of Godsave's analysis and extended and generalized his 
work through a variable property model for the steady burning of a fuel 
droplet in an oxidizing environment. This model included the use of 
integrated forms for the equations of conservation of mass and the 
conservation of energy. In addition, the model included tne 
consideration of the temperature dependence of the fuel thermal 
conductivity and specific heat. While providing a reasonable 
correlation of Godsave ? s experimental results, Goldsmith and Penner 
were also able to formulate explicit expressions for the combustion 
temperature and the flame radius. 

The experiments conducted by Goldsmith [3] to test the theory 
of the Goldsmith and Penner model considered the combustion of 
n-heptane and ethyl alcohol droplets in various ambient oxygen 
concentrations. The theory predicted flame temperatures much higher 
than those found experimentally due to the neglect of dissociation. 

The experimental burning rates compared reasonably well with the 
theory. The experimental range of this study, however, was too limited 
with regard to fuel types, pressure levels, etc. to provide a 
stringent test of the theory. 

Bolt and Boyle [ 18 J obtained the burning rate coefficients for 
droplets burning in a spray but made no measurements of the ambient 
conditions such as temperature and oxygen concentrations. Wood, et al., 
[ 5 J passed droplets through a methane flame and obtained burning rate 
coefficients at 1300°C for a series of fuels under various oxygen 
concentrations, using various drop sizes. They found that the burning 
rate coefficients varied with initial droplet radius and that 


radiative heat transfer had little effect on the results. 
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Faeth and Lazar [ I'J i considered a wide range of bipropellant 
fuel droplets burning in „ . ustion gas environment. This study 
considered an ambient temper i >re range of 1660-2530°K and oxygen mass 
fractions of G-.38. Son* ‘ ; heir findings included the observation 
that all theoretical w u:. - .rcgressively overestimate the burnin c rate 
as the molecular vs, •' the fuel is increased. They attributed 

ud:, failure to * *' ^composition in the region between the droplet 
surface and the ^/^dation zone. They found that the simplified 
diffusion flame theory gave reasonable agreement with experimental 
results . 

Additional studies of low and high pressure convection testing 
may be found in a number of recent reviews [20, 21] on the combustion 
of liquid fuels. 

1. 2.3.2 High Pressure Investigations 

While there has been a large number of combustion studies at 
atmospheric or relatively low pressures, there have been relatively few 
investigations at high pressures. The available high pressure studies 
have shown that certain assumptions of the quasi-steady analysis 
become questionable at very high pressures and particularly near the 
critical point. Williams [22] and Rrzustowski [23] point out that the 
assumptions of negligible finite regression rates of the droplet 
surface and neglecting the transient adjustments of the boundary layer 
around the droplet to changing conditions at the droplet surface 
become suspect as the droplet liquid density approaches that of the 


gas mixture. 
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The experiments of Hall and Diedrichsen [2] considered the 
combustion of suspended droplets in air at pressures up to 20 atm. 

Since this data is presented as total droplet lifetime (which includes 
both heat-up and quasi-steady burning) , the interpretation of the 
results is greatly complicated, 

Brzustowski and Nat.arajan [24] considered the combustion of 
aniline droplets at high pressures in a convective flow field. They 
considered in their investigation a two-film model consisting of a 
spherical inner stagnant film through which diffusion and conduction 
occur and an outer film through which mass transfer occurs. The 
increase in burning rates with pressure whi^h they observed can be 
explained by the effect of natural convection on the combustion 
process. Spalding [25] and Rosner [x6] have considered the 
supercritica 1 combustion of a fuel drop as a problem involving the 
transient diffusion of a dense gas pocket containing pure fuel vapor, 
under conditions of nc convection. Spalding replaces the fuel droplet 
by a point source of fuel vapor, while Rosner considers the fuel vapor 
to be uniformly distributed in a finite-size spherical region. 

Neither of these theories, however, attempts to describe the details of 
the processes occurring in the immediate vicinity of the fuel 
droplet, nor do they define the lower limit of pressure at which they 
are applicable. 

In a theoretical investigation of the transient combustion of 
a liquid fuel droplet which instantaneously reaches its thermodynamic 
critical state on introduction into a supercritical, hot, stagnant 
oxidizing environment, Folymeropoulos and Feskln [26] found that the 
total combustion time decreased with increasing ambient temperatures 
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ambient oxygen concentration. They also found that the flame zone is 
larger than that calculated using models which neglected the 
convection velocity. 

Cnervinsky [27] considered the combustion of liquid droplets 
under supercritical conditions using an analysis in the Von Mises- 
Prandtl plane. He reached essentially the same conclusions as in 
other studies [24, 28, 29] but also found that the burning time is 
dependent upon a temperature term, so that in the case of highly 
exothermic reactions, the burning times are less than those predicted 
by other theories. 

1.2.4 Combustion Under Forced Convection Conditions 
1.2.4. 1 Low Pressure Investigations 

Droplets burn in envelope flames in low velocity gas 
streams whereby the flame envelopes the leading half of the sphere. As 
the velocity increases the boundary layer thickness decreases and 
mass transfer rates increase, but above a critical stream velocity, 
termed the extinction or transition velocity, the flame on the leading 
naif of the sphere is extinguished and a small flame is stabilized in 
the wake of the droplet. [21, 30, 31] Udelson [32] indicates a 
third regime of combustion in which the flame, stabilizes in the 
boundary layer at the sides of the liquid sphere. 

SeveraJ investigations [28, 33, 34] have focussed on the 
variation of the extinction velocity (U^) with droplet diameter (d^) . 
Recent studies [34, 35] have indicated experimentally that the 
extinction velocity is proportional to the square root of the droplet 


diameter . 



There have been few theoretical or experimental investigations 
of the structure of droplet flames Durning under different combustion 
conditions in forced convection. Gollahalli and Brzustowski [35] 
found that the burning rate decreases by a factor of three when the 
envelope flame is transformed into a wake flame at a critical 
approach velocity and that envelope flames radiate much more than 
wake flames. 

Investigations by Sjorg^n [34] and Michael and El-Wakil [36] 
have indicated that in the case of hydrocarbon fuels the formation of 
soot in the wake of the flame may be a significant source of soot in 
spray combustion. Furthermore, this condition is dependent upon 
whether the droplet burns with an envelope or with a wake flame. 

1. 2.4.2 High Fr assure Investigations 

There have been very few theoretical and experimental 
investigations of high pressure, forced convection droplet combustion. 
Several studies [37, 38] have focussed upon the projecting of droplets 
into hot oxidizing environments. Gollahalli [39] studied the flame 
structure in the wakes of small porous spheres burning under forced 
convection conditions. The pressure range was from 1 to 30 atmospheres 
and the ambient air temperature was 300°K. He found that increasing 
the ambient pressure increases the soot concentration and the size of 
soot particles. Also the burning rate decreased by a factor of three 
when the envelope flame was transformed to a wake flame. 

Sami and Ogasawara [40] considered the combustion of fuel 
droplets at pressures up to 16 atmospheres. The ambient air 
temperature was quite low ranging from 25 to 600°C. Their findings 
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included the observations that variations of ambient air temperature 
had little effect on the burning rate and that the burning rate 
decreased markedly when the air velocity exceeded the extinction 
velocity. 


1.2.3 High Pressure Phase Equilibrium Investigations 

Manrique [41] conducted a theoretical study of the vaporization 
of carbon dioxide droplets in a high pressure nitrogen atmosphere. 

His theoretical model considered both steady state and transient 
vaporization cases, which included non-ideal effects associated with 
high pressure mixtures, ambient gas solubility in the liquid phase, 
variation of thermophysical properties in the boundary layer and the 
effects of total pressure on the vapor pressure and enthalpy of 
vapoiization. The ambient temperature range was 375-1600°K and the 
pressure range was 70-120 atmospheres. Some of the conclusions cf 
this work indicate: 

1. A vaporizing droplet can reach and exceed its thermodynamic 
critical temperature, thus becoming a dense mass of vapor, 
at supercritical pressure, by an intrinsically unsteady 
process . 

2. All nonideal effects usually neglected in low pressure 
vaporization models are important in the critical region. 

3. Reasonable estimates of vaporization times over a wide 
range of temperatures and pressures can be obtained with 
the low pressure model provided the effects of ambient gas 
solubility and the effects of total pressure on the vapor 

pressure and enthalpy of vaporization are accounted for 
at high density conditions. 
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The work by Lazar and Faeth [42] and Lazar [43] considered 
bipropellant droplet combustion in the critical pressure range. These 
Let* ts were conducted under zero gravity conditions. The fuels 
considered included n-octar.e and n-decane supported droplets burning 
in air at pressures up to 2000 psia. The theoretical results 
considered the conventional low-pressure treatment of phase 
equilibrium at the droplet surface as well as high pressure models which 
allowed for real gas effects as well as finite ambient gas solubility. 
Some of their findings and major conclusions are as follows: 

1. High pressure corrections and solubility effects should 
be considered in estimating droplet conditions at 

high pressure as well as in estimations of the pressures 
required for supercritical combustion. A good 
approximation for combustion in air can be obtained by 
assuming that the gas at the droplet surface is a 
binary mixture of fuel and nitrogen due to the 
predominance of nitrogen in the system. 

2. The common unity Lewis number assumption will not 
yield useful predictions of droplet conditions at high 
pressures . 

3. Conditions were found where water should condense on 
burning n-paraffin drops. Thi'' occurred for steady 
burning at pressures greater than 1 atm for paraffins 
up to n-pentane. Condensation can also occur for the 
heavier hydrocarbons during heat-up, at sufficiently 


low droplet temperatures. 
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1. 3 Specific Statement of the Problem 

The preceding discussion has indicated that there are certain 
important aspects of high pressure combustion that have rot been 
fully explored. While the results of Manrique, Savery and Borman, 
et al., [13, 41] have demonstrated the importance of high pressure 
effects on the mechanism of evaporation, the combustion of fuel 
droplets was not considered. Other studies [11, 12, 41] also failed 
to consider this aspect of the problem. The investigations of Lazar 
and Faeth [42] provided combustion rate data at high pressures in a 
convection free environment. However, the burning rates obtained were 
unsteady state and the ambient temperature of the air around the 
droplet was low. 

Investigations [39, 40] that considered the effect of forced 
convection on pressurized combustion did not consider realistic 
ambient pressure or ambient temperature levels. This aspect of the 
problem is of great importance when attempting to model an actual 
spray combustion process. 

Therefore, the present investigation was undertaken with the 
following specific objectives: 

1. To compare the theoretical predictions of the various high 
pressure models, when modified to include variable 
properties, ambient gas solubility, and the separate 
determination of the concentrations of the various gas 
phase species with experimental results for hydrocarbon 
and alcohol droplets burning in a high pressure, cold 
air environment under essentially natural convection 
conditions. Comparisons are made between the theoretical 
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and experimental burning rates and liquid surface 
temperatures. Real gas effects are considered and 
emphasis is placed on the prediction of the pressures 
required for the droplets to reach a critical combustion 
condition. 

2. To compare the high pressure theoretical burning rate 
and liquid surface temperature predictions, using the 
above modifications, with experimental results for 
fuel droplets burning at high pressures in hot 
combustion gases of variable oxygen concentrations under 
forced convection. 

3. To determine experimental evaporation rates and liquid 
surface temperatures of fuel droplets, in the absence 
of comhust^n, in hot gases under forced convection and 
to compare the results with the theoretical predictions 
of the various high pressure models. 

The preliminary stage of the investigation consisted of 
developing a high pressure combustion apparatus capable of carrying out 
the experimental objectives. The investigation is in two parts; the 
first aspect of the investigation considered the high-pressure, 
steady-state combustion of liquid fuels in a cold air, natural 
convection environment at pressures up to 100 atmospheres. The 
experimental results were compared with droplet combustion theories 
which both neglected and considered real gas effects. The fuel 
droplets were simulated by porous spheres. Sphere sizes ranged from 
0.63-1.9 cm outside diameter. Fuels considered in this portion of 
the investigation included methanol, ethanol, propanol-1, n-pentane, 
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n-heptane, and n-decane. The ambient gas temperature was maintained 
around 300°K. A discussion of this phase of the investigation is 
presented in Section 4.2. 

The second aspect of the investigation considered the 
combustion and evaporation of fuels in a hot combustion gas environment 
under forced convection. The hot combustion gases were provided by a 
swirl-stabilized burner. The atmospheric tests which provided 
baseline data for the subsequent high pressure vaporization experiment 
employed methanol, ethanol, propanol-1, n-pentane, n-heptane, and 
n-decane as test fuels. The ambient air temperature of these 
experiments ranged from 600°K to 1530°K and over a Reynolds number 
range of 30 to 300. 

For the high pressure tests, the fuels included ethanol and 
n-heptane. The pressure range was from one to 40 atmospheres and 
the ambient gas temperature was maintained around 1150°K for the 
combustion tests and around 1255°K for the pure evaporative tests. 

The Reynolds range for these tests was from 70 to 672. 



CHAPTER II 


THEORETICAL CONSIDERATIONS 

2. 1 Description of the General Model 

The theoretical objective of this investigation was to compare 
the predictions of the various high pressure models with respect to 
their ability to accurately predict bipropellant droplet vaporization 
characteristics at high ambient temperatures and pressures. Major 
empnasis was placed upon the prediction of steady burning rates, 
liquid surface temperatures, and critical burning states under 
variable oxygen concentrations. The present theory is similar in 
many respects to that of References [42 and 43] for high pressure 
droplet combustion. The major point of difference involves the 
different boundary conditions at the liquid surface for the present 
porous sphere combustion case as opposed to steady droplet combustion, 
and a correction to the burning rate for convection effects. These 
considerations are discussed in Sections 2.2.3 and 2.3, respectively. 

The theory may be divided into a gas phase model of the 
combustion process and a phase equilibrium model for conditions at 
the liquid surface. The solution to the theoretical combustion 
problem involves the solution of the conservation equations for 
diffusion of heat and mass. 

The physical model for analyzing quasi-steady bipropellant 
droplet combustion lias be^n developed by Spalding, Goldsmith and 
Penner, et al., [1, 5, and 8). A sketch of this model is shown in 
Figure 1. The present experiment differs from these investigations in 
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that truly steady-state conditions prevail. A major simplification 
to the theory, therefore, is that all transient phenomena (finite 
radius regression rate, time required to attain steady droplet 
temperature, etc.) associated with the combustion process do not have 
to be considered. 

The gas phase analysis involved extending the variable 
property, steady burning solution of Goldsmith and Penner [8] since 
the property variations in this approach are particularly suited to 
the present vaporization conditions. This extension included the 
consideration of the effect of dissolved gas evaporation, separate 
determination of the concentrations of the various gas phase species 
and allowance for variable specific heats of all species. The analysis 
is cast on a molar basis in order to facilitate matching the gas 
phase solution with the phase equilibrium solution at the droplet 
surface. 

The basic analysis, therefore, assumes spherical symmetry and 
neglects convection effects. The analysis of droplet combustion based 
upon the assumption of spherical symmetry with a multiplicative 
correction for convection effects has been demonstrated by many 
investigators [44, 45, 46 and 47] as an approach that gives good 
agreement with experimental results. The treatment of the combustion 
problem under the previously stated assumptions is presented in the 
following sections of this chapter. 

2.2 Gas Phase Model 

2.2,1 Description 


Referring to Figure 1, it is assumed that reaction is confined 
to a spherically symmetric, infinitely thin flame surface where fuel and 
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oxidizer combine in stoichiometric proportions. This assumption 
provides a useful limit when the reaction rates are fast in comparison 
to the diffusion rates, [23] 

The total gas pressure is taken to be constant throughout 
the boundary layer and body forces are neglected. Faeth [48] 
considers this problem and defines the limits under which this 
approximation is valid for droplet combustion. For the conditions of 
the present investigation, consideration of the momentum equation in 
Appendix A.l indicates that the pressure variation across the boundary 
layer is negligible. Thermal diffusion is neglected in the gas phase 
analysis. Williams [44] indicates that thermal diffusion is usually 
negligible in comparison with the ordinary concentration gradient 
diffusion. Therefore, only concentration diffusion is considered in 
the gas phase analysis. 

The influence of compressibility effects on transport 
properties in the boundary layer is neglected and the ideal gas 
equation is employed to compute the molar density in the gas phase. 
Compressibility effects are small except near the droplet surface due 
to the high temperature of most of the gas phase. 

The effect of thermal radiation on the combustion process is 
neglected in the analysis. In the present investigation only 
radiation from the envelope flame is of interest due to the 
construction of the experimental apparatus. The e f fect of radiation 
on the combustion of porous spheres is examined in Appendix A. 2. In 
addition to the above assumptions the concentration dependence of the 
thermal conductivity is neglected and the binary dif fusivities of all 
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species are taken to be the same, although different values of each 
of these properties can be employed inside (Region A, Figure 1) and 
outside (Region b) the flame. The use of equal binary dif fusivi t ies 
in multicomponent mixtures is a common assumption employed in the 
analysis of diffusional processes for burning fuel drop.lets. [44, 49] 

The usual unity Lewis number assumption was not employed in 
the analysis since earlier studies [42] have shown that the value of 
the Lewis number has a strong influence on conditions at the liquid 
surface. The effect of property variations on the combustion process 
is examined in Chapter IV. 

In the present experiment the liquid fuel was pumped to the 
center of the sphere from a storage vessel at atmospheric pressure. 

It is appropriate, therefore, to assume that the liquid entering the 
sphere has a negligible dissolved gas concentration due to the low 
solubility of gases in the test fuels at low pressures. Under this 
assumption, the liquid phase flux of dissolved gas is zero, and the 
fuel is the only component with a finite molar flux inside the flame 
surface . 

The formulation of the problem under the above assumptions 
results in the following unknowns: 

1. Burning rate, n^ 

2. Temperature of the flame, T^ 

3. Liquid surface temperature, T^ 

h. Flame radius, r^ 

5. Species mole fractions at the droplet surface, 

6. Species concentrations in Region B 
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In the following sections of this chapter, and in Appendix B 
fue appropriate transport equation** wre solved to determine these 
parameters for steady droplet turning conditions, 


2.2.2 Conservation Equations 

For a spherically symmetric system (Figure 1) the governing 
equations for the spherical shells defined by Region A and Region 
:>, under the previously states assumptions, for a multicomponent 


system are: 


d'r^V ■ ° l * 1 ' •••’ 

N ^Mass) 

( 2 . 1 ) 

N 

d r 2 V" xt u 2 . 'JT 1 n 

3 — [r ) N.h. - r X-—] = 0 

dr 11 ar 

(Energy) 

( 2 . 2 ) 




dX 



- CD JT + x 1 .I.V 1 ' 1 - 

. . . , N (Species) 

(2.3) 


'lie re is the molar f^ux per unit area of species i. 

In the following the fuel is denoted as component 1 (absent in 
Region B) and the oxidizer, component N (absent in Region A). 

Tiie overall burning rate can be shown to be a constant by 
multiplying Equation (2.1) by 4tt, summing over all species, and 
integrating; this yields 

2 N~1 

4TTr 0 Y ^ = Constant 38 n (2.4) 

i=l 


The mole flux fractions in Region A, are defined as 


N-l 


(2.5) 


Since the fuel is the only component diffusing in Region A, for the 
porous sphere case, and (i-2, N-l) in Region A. 
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Considering the energy equation 


d r 2, . v t 2, dT, r 

— [r (h-h^ - 4Ttr X — ] * 0 


Integrating between r^ and r for r < leads to 


*(h-h t ) - 4 *A a £ 
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ar 


( 2 . 6 ) 


The right hand side of Equation (2.6) may be evaluated through 
consideration of conservation of energy at the droplet interface. 

Under the steady burning assumption .■! of the energy reaching the 
droplet surface goes into the heat of vaporization of the fuel. The 
rate of energy transport is 


, 2, dT 

4TTr X A ^r 


nL, 


(2.7) 


Two cases may be considered theoretically: 

1. For droplet burning, L, Is the heat of vaporization for 
all gas phase components and 


L 



~i L i 


2. For porous spheres, the present case under investigation 
only the fuel is gasified and the quantity must 
include any liquid phase enthalpy rise of the fuel as 
well as the heat of vaporization. 

Therefore, the equation for conservation o' energy for steady 
uurning in Region A becomes: 

h(h-h £ + L) x - 4Ttr 2 X A ^ 


( 2 . 8 ) 
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Employing the definition of given by Equation (2.5) and substituting 
this quantity into Equation (2.3), the conservation of species for 
Region A becomes 

2 dX. 

n[X u -l] - ATTr (CD) A ^i, i-1, ..., N-l (2.9) 

i’he determination of the mole flux fractions in Region B is 
made by considering the combustion process at the flame surface. The 
stoichiometry of the reaction (where W_. denotes the chemical symbol 
of species i) at the flame surface is taken as follows 

N 

W- + Y a.W. (2.10) 

1 & 1 1 


The total molar flux in Region B can be shown to be a constant through 
the use of the continuity equation. This leads to 

N 


n = Z N 

r i=2 

Therefore, for the outer region 

N N 

n = V n, = n, Y a. 
o i 1 i 


i=2 


i=2 


( 2 . 11 ) 


( 2 . 12 ) 


and 
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Writing the conservation of energy across the flame surface 


(2.13) 


(h(h-h £ ) - 47Tr 2 A = {n f £ e i h i - 47rr^X g} (2.14) 


2, dTi 
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The left nand side of the above equation is constant in 
Region A and may be evaluated at the droplet surface. 

For Region B conservation of energy (which must account for 
tne heat of reaction) leads to 
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i^2 1 


(2.15) 


where 


q ji _ h u 




i h i* 


(2.16) 


is the heat of reaction at for gaseous reactants and products. 
The conservation of species in Region B becomes 
N „ dX, 


a i X i " a i } ' 4 * r (CD) B dT 


(2.17) 


2.2.3 Boundary Conditions 

The boundary conditions on Equation (2.8) and (2.9) in 
Region A are as follows: 

r - r r T = T r X ± - X u (i-1, .... N-l) 

r - r f , T - T f , X^^ = 0, x i = x if (i"2 N) 

(2.18) 

In Region B the boundary conditions assume the form: 


r = r f , T = T f , \ m \ f > (W H), \ = 0 

r = r oo» T - T oo* x i - x icc (i=2, .... N) (2.19) 

The steady b lrning equations of :onservation of species and 
energy, Equations (2.8), (2.9), (2.15) and (2.17) and the boundary 
conditions, Equations (2.18) and (2.19), can be integrated upon 
substitution of the relationships for the physical properties 
contained in these expressions. 

For the most general case the f ollowir * relations were 
adopted for the thermal conductivity and the specific heat: 
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c pi * A ± + B^T, i-1, 
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( 2 . 21 ) 
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The quantity 

X. = U/CD). (2.22) 

is only a weak function of temperature and is assumed to have a 
constant average value in each field. 

The detailed solution of the steady burning equations using 
the above property relations is discussed in Appendix B. The results 
for the solution of the unknowns given in Section 2.2.1 are presented 
in the following. 

The most general solution of the equations consider both 
variat\e properties and variable specific heats in the gas phase 
using the property variations given by Equations (2.20), (2.21) and 
(2.22). In the following the fuel is denoted as component 1 and 


oxygen is denoted as component N. 

The gas phase solution in the region inside the flame gives 
the following expression for the fuel mole fraction at the liquid 
surface. There are three branches to the solution: 
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A 1 + Vf + 5 

A i + B i T t - 5 

\ + Vf - 5 | 

A 1 + B 1 h + ?) 


, r > o 

(2.23) 


X U 1 


exp 


2 * a r -i (A i + W -i I A i + B i T f! 1 

* 1 $ I — * — j ] 


4> 2 > o 


(2.24) 


X u * 1 - exp 


2 WW 


(A i + B iV (A i + B i T f } 


<j> =■ 0 (2.25) 


wiiere 


X A ■ WC0) A 



28 


and 

s 2 ■ - * 2 ■ 4 i 2 - 2 B it L i - Vt - r T t 2 > <2 

The mole fractions of the remaining species at the liquid 
surface are given in terms of the composition of the flame 


X iJ£ = X if (1 " X H ) ’ 1=2 ’ N “ 1 (2 

The solution in the region inside the flame yields an 
expression for the burning rate in te.r-'.s r>f properties in the flame 
There are again three branches to the solution as follows: 
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The flame temperature is related to the ambient oxygen 
concentration as follows: 
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where 

Y 2 = - <t' 2 = a' 2 - 2b' (L 1 - Q z - a'T^ - j- T £ 2 ) (2.35) 

An expression for the combustion radius is obtained by 
integrating the energy equation and applying the boundary conditions 
of Equation (2.19). The results are as follows: 


where 
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Another useful approximation is the case when Bj=0, since 
fuel is absent in the outer region and B^ in Equation (2.21) is 
small for the nonfuel species. The most significant difference 
between the solutions for the variable specific heat case and the 
constant specific heat case is the absence of three separate solutions 
for the constant specific heat case. A separate solution is presented 
in Appendix B.2 since the solution for finite does not reduce 
conveniently to this case. The resulting expressions for liquid 
phase mole fractions, burning rate, flame temperature, and flame 
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radius, respectively, are as follows: 


X U " 1 


^ + A L (T f 


V 


1 /A lX A 


Equation (2.25) is unchanged for this case and remains 
X U = X lf (1_ X U>» i=2 > •••» N" 1 
The expression for the burning rate is 
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The expression for the combustion temperature is 
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and the combustion radius is given as 
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The third approximation neglects all temperature dependence 
for all properties. This solution corresponds to the usual low 
pressure, constant property model of droplet combustion. 

The expressions for the liquid phase mole fractions are 
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X 2£ = (1 " X V (2.46) 

The expressions for the burning rate, combustion temperature 


and combvstion radius, respectively, are: 
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The calculations proceeded by guessing a value for the 
liquid surface temperature T^, at a given total pressure, and then 
computing and (the fuel mole fraction at the liquid surface) 

from the phase equilibrium analysis. These values were then employed 
to compute a value for from the gas phase analysis given previously 

and in Appendix B. The final solution was obtained by varying T^ 
until the two values of X^ matched. 


2 . 3 Convection Correction 

Since heat and mass transfer are increased by the effects 
of convection, a low limit for the burninc rate will be obtained if 
the analysis is made for a droplet burning in a still atmosphere 
(stagnant film approximation). In the gas phase model, therefore, 
the effect of convection is treated by the use of a multiplicative 
correction of the burning rate predicted in the absence of convection. 
Since the experiment considers both natural and forced convection, 
the correlations employed are different in each case. 

One aspect of the inves t igatior deals with high-pressure 
droplet combustion under natural convection conditions with air being 
admitted into the combustion chamber at essentially ambient 
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temperature. For these conditions, using the definition of the 
Grashoff number as suggested by Spalding [50] for burning spheres, the 
convection correlation assumes the f orm 

1/1 1 2 

n/n - 1 + 0.221 Pr /J (d g/v) (2.50) 

The properties used in this correlation were taken at the ambient 
conditions of the burning sphere. 

The second portion of the experimental program deals with 
high pressure droplet combustion in a forced convection flow of hot 
combustion gases. Therefore, it is necessary to correct the 
theoretical burning rate constant for the flow of burner gas past the 
test droplet. For steady droplet evaporation or combustion umerous 
investigators [45, 46, and 47] have suggested a convection correlation 
of tne form 


n/n Q - 1 + f (He, Pr, Sc) 

when n Q is the evaporation rate at no-flow. The specific correlation 
that will be used is 

n/n » 1 + 0.278 Re 1/2 Pr 1/3 (1 + 1.237 Re -1 Pr _A/3 ) 

0 

(2.51) 

This expression has the advantage of asymptotically approaching the 
mathematically rigorous results of Fendell, et al., [46] for small 
Peclet number as well as agreeing with the data of Frossling 
[47], Yuge [49], and Allender [50] for 10 < Re < 800. Ambient gas 
properties are used for computing the non-dimensional quantities 
appearing in the above correlation. 



33 


2.4 Phase Equilibrium Analysis 

The determination of steady state condition requires that the 
gas pnase solution given in Appendix B be matched with the boundary 
conditions at the droplet surface. The evaluation of these boundary 
conditions requires a consideration of the thermodynamics of phase 
equilibrium at the vapor-liquid interface of the droplet. Three 
models were employed for computing phase equilibrium at the droplet 
surface, a low pressure model and two phase equilibrium models 
appropriate for use at high pressures. A discussion of these models 
is presented in the following sections of this chapter. 

2.4.1 Low Pressure Theory 

At low pressures it is usually a good approximation to neglect 
solubility effects and take the fuel mole fraction at the droplet 
surface, X^, as the vapor pressure of the pure fuel at the liquid 
temperature divided by the total pressure. Brzustowski [23] has 
snown that this approximation is valid for total pressures up to 1/10 
of the critical pressure of the pure fuel. As the critical pressure 
is approached, high pressure effects cast doubt on the validity of 
the low pressure approximation. 

For the low pressure model, the heat of vaporization, L^, was 

determined by summing the compressed enthalpy change at T q , the heat 

of vaporization at T^, and the ideal gas enthalpy rise between T q and 

T^. In this case the ideal gas enthalpy rise was computed by 

integrating actual specific heat correlations between T and T 0 as 

o x 

opposed to the linearized specific heat correlations employed in 
the gas phase analysis. 
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TUe enthalpy deviations were obtained from the tables of 

Lydersen, et al. , presented in Reference [52], The enthalpy 

deviations were determined between saturated vapor and the ideal gas 

states at T and T fl . 

o Z 


2.4.2 High Pressure Theory 

The two high pressure theories considered solubility and 
other high pressure effects through the use of a modified Redlich- 
Kwong equation of state. In an evaluation of the various methods for 
computing phase behavior, Lazar [43] found that this equation gave 
the most satisfactory agreement with experimental data for binary 
mixtures of a paraffin hydrocarbon with carbon dioxide or nitrogen. 
The Redlich-Kwong equation of state has the form 

o 


P - 


RT 

V-b° 1 0,5 V(V+b°) 


(2.52) 


A vapor phase and a liquid phase are in equilibrium whv,n 
both are at the same temperature and pressure and when the fugacity 
of any conpcnent in the vapor phase is equal to that of the liquid 
phase, i.e., 


f iv " f U* 1-1 N (2,53) 

The fugacity of a component in a mixture is related to the volumetric 

properties of the mixture through the relation 

P 

RT An(f 1 /X 1 P) - f (V ± - RT/P) dP (2.54) 

*'o 

The various component fugacities can be determined by 
substituting the Redlich-Kwong equation of state, using the mixing 
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rules of Prausnitz and Chueh (Appendix C) , into Equation (2.27) and 
integrating. The dimensionless constants ft and ft, required by this 

Si D 

thermodynamic model are obtained by setting the first and second 
isothermal derivatives of pressure with respect to volume equal to 
zero at the critical point. The numerical values of the constants 
then become 0.4278 and 0.0867, respectively. For temperatures far 
removed from the critical region Prausnitz and Chueh [51] suggest 
that it is more appropriate to obtain ft & and ft^ by fitting the 
equation to the volumetric data of the saturated liquid and vapor. 

The binary interaction constants, k^ , are characteristic 
of the i-j interaction for each binary pair present in the system. 

This constant varies between zero and one and increases with increasing 
molecular weight of the hydrocarbon component in systems containing 
hydrocarbons with nitrogen and carbon dioxide. The binary interaction 
parameters, k^ , required by the theory are listed in Reference [42] 
for the paraffins and the combustion product gases. For the alcohols 
the k^j values were taken to be the same as the hydrocarbon 
homomorph of the fuel. The values for the alcohols determined in 
this manner are listed in Table 1. The sensitivity of the final 
calculations to the effect of variation of the interaction parameters 
is discussed in Chapter IV. 

For the high pressure theories, the terms comprising the 
enthalpy rise of vaporization are basically the same as those 
discussed previously for the low pressure theory. The major difference 
between the two theories lies in the computation of the enthalpy 
deviations. As opposed to the method used for the low pressure model. 
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Table 1 


Alcohol Binary Interaction Parameters, 



Substance 

N 2 

co 2 

H 2° 

Methanol 

.10 

.08 

.15 

Ethanol 

.15 

.11 

.20 

Propanol-1 

.20 

.16 

.25 


Note: k . . * k,, and k. . -* 0 

ij ji ii 
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the enthalpy deviations for the high pressure theories were computed 
directly from the Redlich-Kwong equation of state. 

2 . 4 . 2 . 1 Binary Phase Equilibrium Analysis 

For combustion in air the major gaseous species at the 
liquid surface are fuel vapor, nitrogen, carbon dioxide, and water 
vapor. Since nitrogen predominates the non-fuel gases in the system, 
a simplified version of the high pressure theory assumed that this 
system could be presented by a binary mixture of fuel and nitrogen. 

2 . 4 . 2 . 2 Quaternary Phase Equilibrium Analysis 

Previous studies [13,4], 42] have demonstrated the importance 
of gas solubility in the liquid pnase during droplet combustion. 

In the combustion of most bipropellant fuels the principal species 
present in the gas phase include carbon dioxide, water vapor, fuel, 
and nitrogen. Therefore, the most completion version of the high 
pressure phase equilibrium analysis considered this quaternary 


system. 



CHAPTER III 


EXPERIMENTAL APPARATUS 


3 . 1 Introduction 

This chapter is divided into two major parts. The first part 
of the chapter describes the experimental apparatus and procedures 
that were used to obtain fuel burning rates and liquid surface 
temperature measurements for the cold air, natural convection tests. 
The second part of the chapter describes the apparatus and procedures 
used to obtain burning rate and liquid surface temperature data for 
the forced convection tests using 'he high-pressure swirl-stabilized 
burner. 

3.2 Natural Convection Apparatus 

3.2.1 Overall System Descripti on 

A sketch of the overall system is shown in Figure 2. The 
apparatus consists of a high pressure test chamber, gas flow system, 
fuel probe and fuel supply system, and instrumentation to measure gas 
and liquid surface temperatures. 

The test chamber consists of a high pressure cylindrical 
vess ;1 66 cm long with an Internal diameter of 13 cm. The chamber is 
of steel construction and is lined with firebrick. Observation of 
the burning fuel droplets rs made through two quartz windows located 
at the test section. 

The remainder of the systems comprising the experimental 
apparatus is described in detail in th'i following sections. 
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Figure 2 High Pressure Combustion Apparatus (Natural Convection) 
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3.2.2 Gas Flow System 

Air for the combustion process is supplied by a 3000 psia 
compressor. The air flow is accurately controlled by a critical flow 
jeweled orifice system of the type described by Andersen and Friedman 
[52]. The orifices were calibrated using a Precision Scientific wet 
test meter. 

A constant gas composition was maintained around the sphere 
by admitting the air through a multi-holed manifold at the bottom of 
tne combustion chamber. The drift velocity of the air past the 
position of sphere was sufficiently low so that natural convection 
was the predominant flow effect. 

The hot exhaust gases leaving the test chamber were cooled 
in a water cooled concentric tube heat exchanger. Water condensed in 
the heat exchanger was collected in a water trap and periodically 
blown off to a drain. 

The reactor pressure was controlled by a stainless steel 
regulating valve located in the exhaust line. After passing through 
the regulating valve, the reactor gases were exhausted to the 
atmosphere . 

The emergency pressure release system consisted of a rupture 
disc assembly set at 2000 psia. In addition, a one half inch 
stainless steel relief valve was used in conjunction with the rupture 
disc assembly. 
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3.2.3 Fuel Supply System 

The fuel is pumped to the center of the sphere using a 
Whitey precision variable displacement pump. The pump is equipped with 
a 10 mm plunger and has a maximum discharge pressure of 5000 psia and 
a maximum delivery rate of 2200 milliliters per minute. The fuel 
flow rate is measured with a system of graduated burets at the pump 
inlet . 

The porous spheres used in the combustion tests were made of 
alundum. The porous alundum spheres were obtained from the Norton 
Company. Because of irregularities in the manufacturing process, 
it was necessary to grind and sand the rough spheres to a smooth, 
spherical shape. Spheres having diameters of 0.64, 0.95 and 1.90 
cm were employed in the testing. The tolerance on the sphere 
diameters was approximately +0.25 per cent. 

A sketch of a fuel probe is shown in Figure 3. The fuel is 
fed to the center of the sphere through a stainless steel, water 
cooled hypodermic tube and forced radially outward and burned at 
the surface of the sphere. The smallest sphere size was limited 
by the outside diameter of the coolant water tube. The tube 
diameter was 0.20 cm for the 0.64 cm diameter sphere and 0.32 cm for 
the 0.95 and 1.90 cm diameter spheres, respectively. 

Referring to Figure 3, the fuel is introduced through the 
center tube while coolant water passes through the second innermost 
tube and exits through the annulus formed by the second and third 
tubes. The inlet and outlet water coolant temperatures were 
continually monitored by means of thermocouples. For the natural 
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Figure 3 Porous Sphere Probe 
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convection combustion tests the coolant water temperature rise never 
exceeded 2 degrees F. 

3.2.4 Instrumentation 

3. 2. 4.1 Gas Temperature Measurement 

Upstream stagnation air temperatures were measured using 
high pressure Omega chromel-alumel thermocouple probes. The output 
of the gas temperature and liquid surface temperature thermocouples 
was measured with a Leeds and Northrup Model 8686 millivolt 
potentiometer employing an ice bath reference junction. The constant 
temperature reference was provided by a Thermo-Electric Company 
"I-Cell." A series of thermocouples were also located within the 
combustion apparatus to measure the ambient air inlet temperature. 
This value was maintained around 300°K for the natural convection 
tests . 


3. 2.4.2 Liquid Surface Temperature Measurement 

Liquid fuel temperatures at the surface of the sphere were 
made using two 0.0076 cm diameter chromel-alumel thermocouples. 

The thermocouples were mounted on the sphere such that the junctions 
of the thermocouples were flush with the sphere surface. The 
thermocouples were cemented in place with Saureisen high temperature 
cement and located approximately 60 degrees apart along the periphery 
of the sphere. 

3 . 2 . 4 . 3 Dark Field Photography 

A series of dark field photographs of the combustion 


process was taken using a 4 in. x 5 in. Super Graf lex camera fitted 



with a 135 mm Optar Ions. The film usod was Polaroid type 57 
land film with an ASA rating of 3000. The droplet was photographed 
through one of the two observation windows. 

3.2.5 Experimental Technique 

3 . 2 . 3 . 1 Burning Rate Measurement 

Burning rates were determined by measuring the length of 
time required lor a given quantity of the fuel to be consumed. The 
liquid volume measurements were made with a graduated burette 
located at the inlet of the pump. The steady burning rate was 
determined as the flow rate when the surface of the sphere was fully 
wetted and not dripping. Liquid surface temperature measurements 
were also made at this condition. During adjustment to the steady 
burning condition, excess fuel dripping from the sphere was collected 
in a deadended tube at the bottom of the combustion apparatus. The 
drip tube was constructed in such a manner that there was no 
possibility of reignition of the fuel once it entered the tube. 

The liquid fuel was then periodically blown off to a drain. 

Ignition was achieved by bringing a burning match in contact 
with the wetted sphere in the atmospheric baseline tests. The sphere 
was ignited by momentarily placing it in the vicinity of an 
electrically heated nichrome wire in the closed combust! on chamber 
tests . 

3 . 2 . 5 . 2 Operation of the Apparatus 

The preliminary stops in the operation of the apparatus 
involved calibration of the thermocouple output with a millivolt 
potentiometer and focussing the camera for tests where dark field 




photographs were taken. Thor thermocouple-potent iomot or circuit 
was also chocked before runs with an ice bath as reference. The 
fuel pump was then turned on and adjusted so that the sphere was 
fully wetted bui not dripping. The reactor pressure was increased to 
20 psia by admitting air through the manifold. The sphere was then 
brought momentarily in the vicinity of the ignitor wire until 
ignition was achieved. Upon ignition of the fuel probe, the ignitor 
was turned off and the sphere was repositioned bv centering it in the 
observation window at the test section. 

3 . 3 Forced Convection Appa ratus 

3.3.1 Overall Sys t em De s cription 

A sketch of tlie high pressure combustion chamber with the 
burner is shown in Figure 4. The basic reactor vessel is similar 
to the one that was used for the cold gas natural convection tests. 

The reactor consists of a thick-walled cylindrical vessel 66 cm long 
with an inside diameter of 13 cm. The firebrick liner that was 
used for the natural convection tests was not used in the present 
forced convection apparatus. Observations are made through two 
quartz windows located at the droplet test section. 

Other modif icat ions to the basic apparatus included the 
installation of an internal stainless st.'el cooling coil. The 
cooling coil served to reduce radiation from the test chamber walls 
to the test droplet as well as maintaining the chamber walls at 
acceptable temperature levels when under pressure, 

A 10,000 volt spark 'ignition system was installed to ignite 
the burner gases. In addition, a pressurized nitrogen gas system 
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non tod by a 0.7 b kilowatt electric boat or was Installed to romovo 
condensate from the observation windows by directing, a hot jot yf 
nitrogen diluent against thorn. The fuel del Lvorv and measuring system 
was the same as that described for the cold gas tests. 

The basic fuel probe design is the same that was used for 
the natural convection tests. However, because of the hot combustion 
gas environment to which the probes were exposed, the internal 
conf iguration of the probe coolant lines was modified to allow for a 
greater volume flow rate of cooling water. The ext *rnal diameter of the 
probe support was the same as for the natural convection tests. The 
maximum temperature rise of the coolant water was 3 degrees F. 

All of the forced convection tests employed the 0.95 cm 
diameter sphere, both alundum and sintered bronze spheres were used in 
the testing. Subsequent testing revealed no significant differences in 
test results employing the different sphere materials. The bulk of 
the testing was performed employing the alundum porous spheres. 

3.3.1 Swirl-Stabilized Burner 

A sketch of the h gh pressure burner is shown in Figure 5. A 
swirl-stabilized diffusion flame burner is used U\ the apparatus. 

The burner provides accurate composition control of the combustion 
gases, stable operation over a wide range ot flow conditions, and 
the elimination of flashback. In earlier tests with premixed flat 
flame burners, flashback posed serious problems under pressurized 
testing conditions. 

The bur '.or used in the current tests employs a flame supplied 
w*th a mlMure of carbon monoxide and all*. 



Carbon monoxide was chosen 
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Figure S Diffusion Flame Burner fSwlr 1-Stab l l lzed) 
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for a burner fuel since it is quite stable and reduces water vapor 
condensation within the apparatus. The. carbon monoxide used in the 
high pressure burner tests was technical grade and was obtained from the 
Matheson Company, 

The burner is constructed of Norton alundum ceramic tubing and 
type 3 lo stainless steel. The internal diameter of the burner test 
section is 5.08 cm. The burner wall consists of alundum tubing and 
alumina firebrick. Fuel is admitted through a multi-holed injector 
located in the center of the stainless steel burner base. The 
oxidizer is admitted through four tangential injectors located along 
the burner inner wall. A ceramic flow strairhtenor was located in the 
burner passage and was designed such that the velocity and temperature 
profiles of the hot gases were uniform when the test section was 
reached. All gas flows through the burner were metered through a 
critical flow Jeweled orifice system described in Section 3,2.2. 

3.3.3 Instrumentation 

3 . 3 . 3 . 1 Combustion Gas Temperature Measurements 

Ambient gas temperatures at the test section were measured with 
a system of Pt-Pt-10'i Rho fine wire thermocouples . The thermocouples 
were shielded, except for tin last 0.o3 cm, with a 0.079 cm O.D. 
ceramic insulator. The thermocouple measurements when corrected for 
radiation and conduction errors resulted in values that were 
approximately 35 °K higher than the indicated readings. Ambient gas 
temperatures quoted in the investigation are average values for a 


given test condition. 
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To reduce oxidation and catalytic effects the gas temperature 
thermocouples were coated with silicone using a procedure outlined 
by Fristrom and Westenberg [53]. For the ambient gas temperature 
range of the present high pressure tests a silicone coating provided 
adequate protection for the thermocouples. 

As a further check, the gas temperature and composition 
at the droplet location were computed from thermochemical calculations. 
Allowance was made for dissociation and heat loss through the burner 
walls. 


3 . 3 . 3 . 2 Liquid Surface Temperature Measurements 

Liquid surface temperature measurements were made in basically 
the same manner described in Section 3. 1.4.2 for the natural convection 
tests. Because of the hot convective flow environment in which the 
sphere was placed, a third chromel-aluroel thermocouple was located at 
the top of the sphere and flush with the surface to see if there were 
any significant variations in liquid surface temperature around the 
periphery f the sphere. In addition to this technique the sphere 
was rotated at intervals up to 360 degrees from its* original 
position to provide a further check. Only minor variations were 
noted in the liquid surface temperatures. Consequently, only the 
two thermocouples located on the lower half of the sphere surface 
were used to obtain test results. 

The two chrome 1-alumel thermocouples were located closer to 
the bottom stagnation point of the sphere (approximately 30 degrees 
apart) since preliminary testing with the burner revealed a tendency 
for the sphere to dry off at this point at high combustion gas flows 


\ 
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and at high ambient temperatures Locating the thermocouples in 
tnis manner allowed this condition to be easily detected and 
corrected. 

3 . 3 . 3 . 3 Operation of the Apparatus 

Preliminary steps in the operating procedures included 
calibrating and checking the liquid surface temperature and combustion 
gas temperature thermocouples. The combustion chamber and burner 
were then purged with nitrogen. The burner ignition sequence 
consisted of activating the spark ignitor and then turning on the carbon 
monoxide-air mixture to the burner. 

Following ignition the ignitor uas turned off and the gas 
flows carefully adjusted. The burner combustion gas temperature and 
the prob r and cooling coil water temperatures were continually 
monitored throughout the test. The fuel probe was ignited directly by 
the hot combustion gases. The recording of data followed essentially 
the same procedure as that used for the free convection tests. 











CHAPTER IV 


THEORETICAL AND EXPERIMENTAL RESULTS 

4 . 1 Introdu c tion 

One of the major objectives of the present investigation was 
to determine bipropellant droplet vaporization characteristics under 
essentially combustion chamber conditions. Since convection has a 
dominant effect on the burning rate of fuel droplets, the experiment 
focussed on the combustion of fuel droplets under natural and forced 
convection conditions. The tests extended over a wide range of 
ambient pressures , ambient gas temperatures, Reynolds numbers, and 
ambient oxygen concentrations. Also, a wide range of fuels and 
droplet sizes was considered. This chapter is divided into two 
parts. The first section discusses the results o r the natural 
convection tests, while the second section considers the results of 
forced convection tests in a combustion chamber environment. 

4.2 Natural Convection Tests 

4.2,1 Fuels and Range of the Test s 

The objective of this part cf the investigation was to study 
high pressure droplet combustion for a vatiety of fuels in a cold gas, 
natural convection environment. The experimental determination of 
both burning rates and liquid surface temperatures was considered In 
tt\is effort. The experimental results were compared with droplet 
combustion theories which both neglected and considered real gas 
effects. The fuels used in the test consisted of three alcohols and 
three n-paraffins. The alochol3 included methanol, ethanol, propanol-1 
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the n-paraf f ins included n-pentane, n-heptane, and n-decane. The 
pressure range varied from 1 to 51 atmospheres for the n- paraffin 
fuels and from 1 to 78 atmospheres for the alochols. korous 
alundum spheres were used to simulate the fuel dropJets. The sphere 
sizes included 0.63 cm, 0.95 cm, and 1.9 cm di-aiie^cr spheres. 

In addition to the determination of burning rates and liquid 
surface temperatures, a series of dark field photographs were taken. 
Methanol and ethanol were the fuels used in these tests. The sphere 
sizes employed were 0.63 cm, 0.95 cm, and 1.9 cm, in the pressure range 
from 1 to 35 atmospheres. 

The experimental burning rate measurements and liquid surface 
temperature measurements reported in the following section of this 
chapter for the natural convection tests represent the average 
values for a series of separate runs. Typically, the values from the 
individual tests were within 2% of the average. 

4.2.2 Observations 

Envelope flames were observed during the natural convection 
tests with alcohols. The flames were blue over the forward portion 
of the sphere surface, turning to a yellowish flame over the rear half 
and in the wake. A photograph of a burning methanol droplet at 20 psia 
using the 0.95 cm sphere is shown in Figure 6. The same droplet 
burning at an ambient pressure of 310 psia is shown in Figure 7. At 
the higher pressure the envelope flame is very near the droplet surface 
and a luminous wake extands several sphere diameters above the probe. 

The combustion flames of the n-paraffin fuels were very 
luminous. Sooting was a serious problem with all of the n-paraffin 
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fuels. Methanol gave no indication of soot, and ethanol only began 
sooting at very high pressures (around 1100 psia) . Propanol-1, 
on the other hand, exhibited sooting problems at much lower pressures 
(around 40 atmospheres). 

Sooting was first observed around 7 atmospheres for decane 
and heptane and around 6 atmospheres for n-pentane. The upper limit 
for data acquisition due to sooting was approximately 52 atmospheres 
for n-pentane and n-heptane and around 32 atmospheres for n-decane, 
although the test renge extended to 55 atmospheres. 

4.2.3 Burning Rates 

All the natural convection experimental results were 
obtained for combustion in air. The ambient air temperature and the 
fuel inlet temperature were both 300 K for these tests. The 0.95 cm 
diameter alundum sphere was used for the bulk of the burning rate 
measurements . 

Figures 8 and 9 illustrate the burning rate results for the 
alcohols and paraffins, respectively. The theoretical results shown 
in Figures 8 and 9 were calculated ^ing the variable property- 
variable specific heat gas phase ai~ sis. The quaternary version is 
illustrated for the high-pressure theory. 

The theoretical curves are terminated at high pressures 
when the critical burning condition is reached. For the low-pressure 
theory, critical burning was assumed to occur when the liquid surface 
temperature was equal to the critical temperature of the fuel. 

Critical burning for the high pressure theory formally occurs when the 
liquid surface reaches its critical mixing point for the conditions of 
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Figure 8 Experimental Burning Rates for Alcohols 


BURNING RATE (g /sec.) 
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PRESSURE (atm) 

Figure 9 Experimental Burning Rates for N-Paraffins 
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the combustion process. The burring race predictions of the two 
theories are almost identical, although the high pressure theory 
generally predicts a higher pressure for critical combustion. 

The experimental results for methanol and ethanol (Figure 8) 
were terminated at high pressures due to difficulties in determining 
th*± burning rate. At pressures on the order of 80-100 atm, for these 
fuels, the flame zone would tend to move away from the sphere with 
increased fuel flow rates and clear evidence cf fuel dripping could not 
be ob earned. For both of these fuels vapor jets were observed rather 
than liquid drops, although methanol presented this problem at 
slightly higher pressures. This behavior probably indicates the 
onset of critical burning for these fuels, but the pressure at which 
this condition occurred could not be defined very precisely. 

The burning rates for the remaining fuels in Figures 8 and 
9 are terminated at high pressures due to the formation of soot. Ir 
these cases, carbon spots would form and grow on the surface of the 
sphere causing the test to be terminated at elevated pressures. The 
upper pressure limit for testing decreased with increasing molecular 
weight for the paraffin fuels tested. 

The absolute agreement between the theoretical and experimental 
burning rates in Figures 8 and 9 is comparable to results obtained by 
Faeth and Lazar [19] in low pressure tests. In particular, the theory 
gives a reaonsably good indication of the rate of increase of the 
burning rate with increasing pressure. 

The effect of varying sphere size is examined in Figure 10. 

The sphere sizes employed ranged from 0.63 cm-1.9 cm. For this plot, 



nT,b,/4Ilr.X t )/(n / n 
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PRESSURE (otm) 

Figure 10 Burning Rates for Variable Sphere Size 
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the dimensionless burning rate, normalized by the convection correction 
is employed for the ordinate so that data for various sphere sizes 
should fall on a single curve. 

The results shown in Figure 10 indicate that the normalized 
burning rate (which corresponds to the no-convection burning rate of 
the theory) is almost a constant up to the critical burning condition 
for the present porous sphere experiments. This is due to the fact 
that the no-convection burning rate is largely dependent upon the 
total enthalpy rise of vaporization, which does not change to a great 
degree with increasing pressure for porous spheres. For porous sphere 
combustion, the reduced heat of vaporization near the critical point 
is compensated by increases in the enthalpy rise required to bring the 
fuel from the inlet to the surface temperature. 

The fact that the normalized burning rate is relatively 

constant indicates that the increase in burning rate with increasing 

pressure in Figures 8 and 9 is largely due to convection effects. The 

present experimental results represent a reasonably good test of the 

burning rate correction for natural convection, since the Grashoff 

number, based upon the Spalding [55] definition used in Equation (2.50), 

4 8 

varies in the range 10 -10 . 

4.2.4 Liquid Surface Temperatures 

The liquid surface temperature results for the six fuels are 
illustrated in Figures 11 and 12. The boiling point curves and the 
surface temperature predictions of both the low-pressure and the 
quaternary high-pressure theories are shown on the figures along with 


the data. 
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The difference between the two theories is inure obvious 
with regard to surface temperatures, than was the case for the burning 
rates, with the high pressure theory predicting the lowest surface 
temperature at a given pressure. 

It is seen in Figure 11 that the data for ethanol and 
propanel-1 agrees reasonably well with the high pressure theory at 
high pressures. For methanol, however, the low pressure theory 
gives the best estimation of the data over the entire test range. 

The poorer high-pressure theoretical results for methanol could be due 
to the large quantities of jwater vapor in the combustion products of 
this fuel. Water is difficult to model precisely in the high-pressure 
phase equilibrium analysis, and material-" with high water vapor 
concentrations in the products have generally shown poorer agreement 
with the high-pressure theory in the past. [56] 

The experimental liquid temperature data for the paraffins 
shown in Figure 12, could not be extended lO sufficiently high 
pressures to provide an adequate test of the high pressure theory due 
to the formation of soot. Over the available experimental range, the 
low pressure theory appears to ne adequate for these materials. 

4.2.5 Phase Equilibrium Results 

Figure 13 illustrates computed gas and liquid phase 
compositions, at the liquid surface, for propanol-1 and n-heptane. 
These results pe v .ain to porous sphere combustion in air, with a fuel 
inlet and ambient a^r temperature of 300°K. The gas phase composition 
remains relatively constant as the total pressure is increased for 
both fuels. In contrast, the liquid phase concentration of dissolved 
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PRESSURE (atm) 


Figure 13 Predicted Liquid Surface Compositions for Porous-Sphere 

Combustion in Air. with Fuel Inlet and Ambient Temperature 
of 300*K 
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gas increases significantly with increasing pressure. The critical 
mixing point of the surface (the critical combustion condition) is 
indicated by the equality of the liquid and gcc phase composition 
at this state. The dissolved gas concentration becomes quite large 
near the critical combustion condition for the present test 
conditions, reaching values as high as 60% for n-decane. 

At pressures higher than the critical combustion condition, 
the process is similar to the porous sphere combustion of a gas. In 
this regime, no liquid surface would be observed and a range of fuel 
flow rates (subject to blow-off and quenching limits) could be 
accommodated by the sphere at a given pressure, as opposed to th< single 
fuel flow rate possible for liquid fuel combustion at pressures below 
the critical combustion condition. 

4.2.6 Discussion 

The preceding theoretical results were obtained with the 
quaternary phr.se equilibrium model. The simplified binary model gave 
essentially the seme results with regard to burning rates and liquid 
surface temperatures. In contrast to high pressure droplet combustion, 
however, there were significant differences between f he critical porous 
sphere combustion pressures predicted by the two high-pressure 
theories. The critical combustion conditions for all three theories 
are compared with pure fuel critical properties in Table 2. In 
agreement with the experimental findings, both the low pressure and 
high pressure quaternary theories predict critical burning pressures 
on the order of ?00 atm for methanol and ethanol. The theoretical 
indication that oth&nol should experience critical burning at 
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pressures somewhat below methanol is also in qualitative agreement with 
the fact that experimental difficulties in determining burning rates 
were encountered at somewhat lower pressures for ethanol, c.f. Figure 8. 

All the previous theoretical results were obtained with the 
variable property — variable specific heat gas phase analysis, using 
the properties listed in Table 3. The use of the variable property- 
constant specific heat and constant property gas phase analysis 
gave essentially the same results, when the respective constant 
properties in each of these cases were evaluated at average conditions 
in each region. The effect of parametric variations of the and 
the gas phase properties listed in Table 1 was also examined. The 
value of had the greatest influence on the prediction of liquid 
surface temperatures and critical burning conditions. Quantitatively, 
the effect of variations of this parameter was similar to that 
encountered in earlier studies of high pressure combustion. [56, 57] 
Variations in the predicted burning rates were almost in direct 
proportion to variations in the value of A^, and were relatively 
insensitive to changes in 


4 . 3 Forced Convection Tests 

4.3.1 Fuels and Range of the Tests 

The objective of this part of the invest? ga‘ ion was the 
determination of high pressure droplet burning and evaporation rates 
and liquid s rface temperatures in a simulated combustion chamber 
environment. The experiment considered various gas flow velocities 
past the test droplet and various pressure levels to provide an 
evaluation of both the effects of forced convection and high pressure 
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*At 1000°K., X from Equation (2) at any other temperature. 
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phenomena on droplet combustion. The experimental tests focussed 

upon droplet combustion and evaporation at atmospheric and high ambient 

pressures. 

The tests at atmospheric pressure provided baseline data for 
the high pressure tests. Fuels considered in this phase of the 
investigation included methanol, ethanol, propanol-1, n-pentane, 
n-heptane, and n-decane. The ambient gas temperature range was from 
600 to 1530°K and the ambient oxygen molar concentrations included 
13%, 9.5% and pure evaporation. The Reynolds number range, based 
upon approach conditions, varied from 30 to 300. These baseline 
tests employed the swirl-stabilized burner described in Chapter III 
to provide the hot combustion gas environment around the test droplets. 
The 0.95 cm sphere was employed in the tests, using both alundum 
and sintered bronze spheres. 

For the high pressure droplet tests, the fuels included 
ethanol and n-heptane. Three ambient molar oxygen concentrations 
were employed. These included molar concentrations of zero, 9.5 and 
13 percent, respectively. The two combustion cases considered ambient 
oxygen concentrations of 9.5 and 13 percent and the ambient gas 
temperature was maintained around 1150°K. For the purely evaporative 
case the ambient gas temperature was maintained around 1260°K. 

The Reynolds number range for the elevated pressure tests 
was from 70 to 672 and the pressure range was from one to 40 atmospheres. 
The tests employed a porous sphere having a 0.95 cm diameter. Eoth 
alundum and sintered bronze spheres were, used, although the bulk of 
the measurements were made using the alundum spheres. 
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The experimental burning rate measurements and liquid surface 
temperature measurements reported in the remainder of this chapter for 
the atmospheric baseline and high pressure forced convection tests 
represent the average values for a series of separate runs performed 
at the same ambient gas temperatures and ambient pressures. The 
values obtained from the individual tests were within 3% of the 
average. 

4.3.2 Observations 

A pale blue flame was observed around the leading face of 
the sphere and close to the sphere surface for the combustion of the 
alcohols. The tail of the flame was yellowish and extended several 
diameters above the sphere. These characteristics are typical of an 
envelope flame. The combustion flames of the n-paraffins were almost 
completely yellow, and the flame of n-decane, in particular, was 
quite smoky. 

For a finite ambient oxygen concentration, a well-defined 
luminous wake was present around the tect droplet. With decreasing 
ambient oxygen concentration, the flame zone would move away from 
the droplet and the intensity of the luminosity would decrease. For 
negligible ambient oxygen concentration (evaporative case), a 
diffuse luminous wake was present at high ambient gas temperatures 
even though there was no flame resulting from exothermic reaction. 

Tnis wake was probably due to radiation from hot carbon particles 
formed by decomposition of the fuel in the burner gases. 

In agreement with the observations of other investigators [19], 
methanol did not exhibit a luminous wake at negligible ambient oxygen 
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concentrations and high ambient gas temperatures. This behavior may 
be attributed to the high resistance of methanol to the formation 
of carbon particles upon decomposition. The behavior of methanol was 
similar to the other fuels under oxidation conditions, however, 
although the flame luminosity was less intense. 

At higher pressures the luminosity of the combustion flames 
for both ethanol and n-heptane increased dramatically with increasing 
ambient oxygen concentration. The rate of soot formation intensified 
with decreasing ambient oxygen concentratioi as the pressure was 
increased for the combustion cases employing n-heptane as a fuel. 

In the open atmosphere baseline tests, side flames were 
observed at very high Reynolds numbers. These flames are characterized 
by the flame stabilizing along the sides of the sphere with the 
forward portion of the sphere extinguished. This phenomenon was 
not observed in the pressurized tests. No attempt was made to 
determine precisely when side flames occurred, since the determination 
of blow-off conditions was not an experimental objective. All 
measurements were made when the droplet burned with an envelope 
flame. There was no evidence of strictly wake flames, in which the 
flame resides totally in the downstream wake of the sphere, over the 
Reyrolds number, ambient pressure and ambient gas temperature range 
of the tests. 

4.3.3 Burning and Evaporation Rates 
4 * 3. 3. 1 Atmospheric Pressure Baseline Test Results 

Table 4 lists the computed properties of the ambient gas 
at the droplet test section for the most widely run test conditions. 
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The oxygen mole fraction in the ambient gases quoted in this table 
is an effective mole fraction based upon the concentrations of 
possible oxidizing species (U^, 0, NO) since the minor species 
(0, NO, etc.) are present only in very small concentrations at these 
gas temperatures. 


Table 4 


Properties of the Ambient Gas for the Test Conditions 



X 02°° 

Too 

°K 

CO 

co 2 

N 2 


0.130 

1145 

0.000 

.132 

.738 

0.095 

1145 

0.000 

.190 

.715 

0.000 

1255 

0.000 

.348 

.652 


Plots of the experimental and theoretical burning rates versus 
ambient gas temperature at atmospheric pressure are shown for the 
alcohols in Figure 14 , while Figure 15 gives similar results for the 
n-paraffins. The experimental results v/ere corrected by the use of 
the multiplicative correction for forced convection given by Equation 
(2.51). The theoretical predictions were made using the variable 
property-variable specific beat low pressure model of Lazar and Faeth 
[42]. The theoretical burning rates tended to be lower than the 
measured rates for all fuels. The theory does predict the correct 
trends in vaporization rates with ambient gas temperature as may be 
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Figure 14 Atmospheric Pressure Burning Rates of the Alcohols, Ambient 
Oxygen Concentration 9.5% (Molar) 
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Figure 15 Atmospheric Pressure Burning Rates of the N-Paraffins, Ambient 
Oxygen Concentration 9.5% (Molar) 
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noted from the slopes of the curves. These tests results are in 
contrast to the burning rate measurements obtained for the high 
pressure, natural convection, cold gas tests discussed in the previous 
section. For these tests the theoretical predictions were found to 
be, on the average, higher than the observed burning rate. The 
experimental results did, however, confirm the theoretical prediction 
that ambient gas temperature has little effect on the burning rate. 

Evaporation rates with zero oxygen concentration in the 
ambient gases are shown in Figures 16 and 17 for the different fuels. 
The theoretical predictions were lower than the experimental values 
over the entire test range. The evaporation rates followed a fairly 
linear increase with ambient gas temperature. 

The evaporative case for methanol in Figure 18 is shown over 
a wider temperature range than the earlier results. From this plot 
it can be seen that the correct trend in evaporation rates is observed 
between theory and experiment and that there is fair agreement at the 
lower ambient gas temperatures. A consideration of the evaporation 
rates in Figures 16 and 17 reveals that an extrapolation of the 
experimental and theoretical results in these figures lead to the same 
general conclusion. 

The behavior of the theoretical evaporation curves at various 
ambient gas temperatures is examined in Figure 19 for etnauol and 
n-heptane. The results are for vaporization rates from a 0.95 cm 
sphere at one atmosphere pressure. The results show that with 
decreasing ambient gas temperatures, the vaporization rates are quite 
sensitive to ambient oxygen concentration. 
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Figure 17 Atmospheric Pressure Evtpora’j >n Rates for N-Paraffins 
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Figure 18 Methanol Evaporation Rates 
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Liquid surface temperatures as a function of gas temperatures 
for ethanol and n-pentane are shown in Figure 20. The low pressure 
theory tended to predict lower droplet surface temperatures than the 
.^asured results. This was observed for both the combustion and 
evaporative cases, although the difference was more pronounced with 
respect to the evaporative case. 

The effect of ambient oxygen concentration on the vaporization 
rate at a fixed ambient gas temperature of 1145°K is examined in 
Figure 21 for ethanol and n-heptane. There is fairly good agreement 
between theory and experiment at higher ambient oxygen concentrations. 
The poorer agreement between experiment and theory at low ambient 
oxygen concentrations was observed by other investigators [19] under 
similar test conditions. 

4. 3. 3. 2 High Pressure Test Results 

Plots of the experimental burning rate versus pressure 
at various ambient oxygen concentrations are shown in Figures 22 and 
23 for ethanol and n-heptane. The theoretical predictions of the 
variable property-variable specific heat versions of the low pressure 
and high pressure quaternary theories are shown along with ♦‘he data. 

For these tests the ambient gas temperature was maintained around 
1145°K (1600°F) for the two combustion cases and 1255°K (1800°F) for 
the evaporative case. The sphere size employed for these tests was 
0.93 cm in diameter. 

The theoretical curves are terminated at high pressures 
when the critical burning condition is reached (denoted by an asterisk) . 
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for Ethanol and N-Pentaue 
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The burning rate predictions of the two theories are in * 

reasonable agreement with the experimental results at the higher 
ambient gas temperatures considered in these tests. These findings 
are consistent with the results of the natural convection, cold gas 
test results. The greatest difference between the two theories lies 
in the predicted critical burning condition. The low pressure theory 
predicts a significantly lower pressure for critical combustion, at 
a given ambient oxygen concentration, than the high pressure theory. 

The difference between the estimations of critical combustion pressures 
is greater in the case of n-heptane than for ethanol. Both theories 
predict an increase in the ambient pressure required for critical 

burning as the ambient oxygen concentration is reduced. 

The experimental curves for n-heptane were termined due to 
the formation of soot for the two combustion cases. The combustion 
flame for the .095 molar oxygen concentration case was particularly 
smoky at pressures in the vicinity of 20 atmospheres. The experimental 
curve for the evaporative case was terminated due to experimental ; 

difficulties associated with maintaining a fully wetted sphere surface. j 

! 

At high pressures and high Reynolds numbers, there was a tendency \ 

for the sphere to dry off at the bottom stagnation point and then to j 

progressively overheat. j 

/ 

The experimental curves for ethanol were terminated mainly j 

due to difficulties in maintaining a fully wetted sphere at high j 

pressures. There was no evidence of sooting over the ambient oxygen j 

and pressure range considered in the tests. In both the ethanol and \ 

n-heptane combustion cases, the luminosity of the flame increased with i 

* 

increasing ambient oxygen concentration. 

i 

\ 

— - - WTT STm'S?'****** ft 
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The vaporization rates shown in Figures 22 and 23 reveal a 
similar trend in burning rates with pressure noted for the normalized 
burning rates (Figure 10) for the natural convection tests, ihe 
present test results represent an adequate test of the burning 
rate correction for forced convection (Equation (2.51)) since the 
test results (which were taken over a Reynolds number range of 75 to 
670) cover a significant portion of the Reynolds number range for 
which this correlation is applicable. 

The liquid surface temperature results for ethanol and 
n-heptane are illustrated in Figures 24 and 25. The boiling point 
curve and the surface temperature predictions of both the low-pressure 
and quaternary high pressure theories are shown on the figures 
along with the data. The difference between the two theories is 
more obvious with regard to the prediction of surface temperatures. 

For the liquid temperature results, the same general trends 
that were observed for the cold gas, natural convection tests can 
be seen for the high ambient gas temperature, forced convection tests. 
The n-heptane measurements tend to follow the low-pressure predictions, 
whereas the ethanol experimental results did not extend to high 
enough pressures to provide an adequate test of the theory. Over the 
experimental range shown, the low pressure theory appears to be 
adequate . 

The theoretical predictions of the binary high pressure 
theory are examined in Table 5. This table considers the critical 
combustion pressures for ethanol and n-heptane at an ambient gas 
temperature of 1145°K. While both the high pressure binary and 
quaternary theories predict approximately the same liquid surface 






Figure 25 N-Heptane Liquid Surface Temperatures 
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temperatures, there is a significant difference in the predicted 
critical combustion pressures. The experimental tests ii. this 
investigation considered a four-fold variation in ambient gas 
temperatures (300°K to 1255°K) at high pressures. Over this ambient 
temperature and pressure range, there is a consistent tendency 
between the predictions of the two high pressure versions, with the 
binary theory generally predicting a much higher pressure for 
critical burning conditions. 


Table 5 


High Pressure Theoretical Predictions for Critical Burning Conditions 


Fuel 

Ethanol 


N-Heptane 

Theory 

x 0 2 - * 095 
t r 

P R 

T 

R 

P R 

binary 

.949 

2.22 

.943 

4.22 

Quaternary 

.946 

1.83 

.936 

3.26 

binary 

x o 2 * * 13 

.957 

2.14 

.945 

3.89 

Quaternary 

.953 

1.67 

.939 

3.04 


The predicted reduced pressures for critical burning at 
various am' lent oxygen concentrations and ambient gas temperatures using 
the low pressure variable property-', .riable specific heat model are 
shown in Table 6 for ethanol and n-huptane. The pressure required for 
critical burning decreases with increasing ambient oxygen concentration, 





91 


Table 6 


Predicted Reduced Pressures for Critical burning (L.P. Theory) 
at Various Ambient Gas Temperatures 


Ambient Gas Temperature (°K.) 

Fuel X 1145 1255 1365 1480 1645 1920 


Ethanol .210 

.130 
.095 
.005 

.210 
.130 
.095 
.005 


1.37 

1.35 

1.35 

1.49 

1.48 

1.46 

1.51 

1.51 

1.51 

1.97 

1.97 

1.97 

1.78 

1.74 

1.74 

1.85 

1.85 

1.85 

1.85 

1.85 

1.85 

1.85 

1.85 

1.85 


1.33 

1.32 

1.30 

1.44 

1.41 

1.38 

1.51 

1.51 

1.51 

1.97 

1.97 

1.97 

1.70 

1.70 

1.67 

1.81 

1.81 

1.74 

1.85 

1.85 

1.81 

1.85 

1.85 

1.85 


n-Heptane 
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altuough the difference is not appreciable considering the ambient gas 
Lemperalure range involved. 

4 . 3 . 3 . 3 Phase Equilibrium Results 

Figures 26 and 27 illustrate computed gas an liquid phase 
compositions at the liquid surface for n-heptane and ethanol. The 
results in Figure 26 are for porous sphere combustion in air with a 
fuel inlet temperature of 30G°K and an ambient gas temperature of 
1145 3 K. There are no significant differences between the results for 
the two ambient oxygen concentrations shov:n in this figure. 

Trie binary th retical predictions for ethanol are shown in 
Figure 2‘ . The predicted gas phase mole fractions for the fuel and 
nitrogen were more sensitive to ambient oxygen concentration than 
was the case for the quaternary theoretical results. 

4. 3.3.4 Discussion 

The low pressure and high pressure theories gave essentially 
the same predictions of burning rates at high pressures in agreement 
with the cold gas, natural convection tests. The convection 
correction for forced convection gave adequate results for the 
Reynolds number range 70-672 over which the tests extended, although 
there is poorer agreement betweer theory and experiment for the 
evaporative case. Similar findings were reported in Reference [19] 
for low pressure tests. Despite the discrepancy between the theoretical 
and experimental results, both predict the same trends in evaporation 
rates with respect to ambient gas temperature and pressures * 

Essentially the same general conclusions that were reached 
for the natural convection tests are applicable to the forced 
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Figure Predicted Liquid Surfec* Compositions for Porous-Sphere 
Combustion jith s Fuel Inlet Temperetare >f 300*K end en 
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Figure 27 Predicted Liquid Surface Compositions for Porous-Sphere 
Combustion with a Fuel Inlet Temperature of 300°K and an 
Ambient Gas Temperature of 1145 °K. (Binary Theory) 
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convection tests. Even at the higher ambient gas temperatures , the 
low and high pressure models predicted the same trends that were 
noted for the cold gas tests. These included the lower liquid surface 
temperature predictions of the quaternary model as compared to those 
of the low pressure model. 

Another similarity between the theoretical results of the 
natural and forced convection tests occurred in the predictions 
obtained using the variable property constant specific heat gas 
phase analysis. For the cold gas tests both the variable property- 
variable specific heat analyses gave essentially the same results as 
long as the respective constant properties were evaluated at average 
conditions inside and outside the flame* For the high ambient gas 
temperature tests both the quaternary and binary high pressure models 
gave essentially the same results for the variable and constant 
specific heat cases, although the binary model consistently predicted 
much higher pressures for critical burning than the quaternary theory. 

On the other hand, the low pressure model gave identical 
predictions tor both the variable specific heat and the constant 
specific heat cases at all ambient oxygen concentrations. Another 
interesting observation concerning the low pressure model is that 
over an ambient gas temperature range of 300°K to 1925 °K, the predicted 
critical combustion pressure for ethanol at 21 per cent ambient 
oxygen molar concentration varied from only 88 to 82 atmospheres, 
respectively. The high pressure theories predicted a much wider 
range of pressures for critical combustion conditions at a given 
ambient oxygen concentration for a similar temperature range. 
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The sensitivity of the calculated results to errors in the 

properties was examined by parametrically varying x A (+20%) and the 

binary interaction parameter k _ (+20%) . The effect of varying the 

modified Lewis number (x^) for ethanol and n-heptane during combustion 

(X Q = .095) at an ambient gas temperature of 1145 °K is examined in 
2 

Figure 2£. The theoretical results were obtained using the variable- 
property-variable specific heat low pressure model and the high- 
pressure quaternary model. In agreement with the cold gas test 
results the variations of x A (+20%) had the greatest effect on the 
predicted critical burning conditions. The low pressure predictions 
are virtually the same for n-heptane but the predicted critical 
pressures for ethanol vary an average of 15 per cent for the conditions 
shown. 


The burning rate predictions obtained for v * low and high 
pressure models were essentially the same as those predicted for the 
normal case. 

The variation of the binary interaction parameter k^ (+20%) 

had very little effect on the predicted critical burning liquid 

surface temperature. The temperature variation and reduced pressures 

for critical burning (X - .095) for ethanol and n-heptane is 

2 

examined in Table 7. The percentage variation is less than 2% 
for the predicted liquid surface temperatures and less than 6% 
for the predicted critical burning pressures. The results shown 
were obtained using the high-pressure variable property-variable 
specific heat quaternary theory. A similar trend was noted using 
variable oxygen concentrations and employing the high pressure binary 
theory. 
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Figure 28 Effect of Variation of the x A Parameter on Droplet combustion 
(X - .095) 
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Table 7 


Effect of Variation of the k-y Parameter on Predicted Critical 
Burning Conditions (X^ = .095) 


Fuel 


Ethanol 


n-Heptane 

Variation 

t r 

P R 

t r 

P R 

Normal 

.945 

1.825 

.936 

3.26 

+20% k. , 
ij 

.941 

1.944 

.934 

3.33 

-20% k ii 

.949 

1.714 

.925 

3.02 





CHAPTER V 


SUMMARY AND CONCLUSIONS 


5.1 Summary 

The overall objective of *“he present investigation was to 
determine the characteristics of bipropellant droplet vaporization 
at high pressures under natural and forced convection conditions. 

The tests considered both cold gas and high ambient gas temperatures. 
The high pressure forced convection tests were conducted under 
simulated combustion chamber conditions. Porous spheres were used to 
simulate the fuel droplets in order to determine steady droplet 
burning rates and liquid surface temperatures. The tests considered 
various ambient oxygen concentrations and sphere sizes. Major 
emphasis was placed upon studying the effects of natural and forced 
convection on the combustion process. 

The theoretical model that was used in the investigation was 
developed by extending the variable property, steady burning theory of 
Goldsmith and Penner [17]. The extensions included the effects of 
dissolved gas evaporation, allowance for variable specific heats of 
all species and the separate determination of the concentrations of 
the various gas phase species. The effect of convection was treated 
in the analysis through the use of a multiplicative correction of 
the theoretical burning rate computed for the case of no convection. 

In the phase equilibrium analysis ambient gas solubility 
and high pressure effects were considered in the determination of 
conditions at the droplet surface. The calculations for these 
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conditions employed the modified Redlich and Kwong equation of state 
along with suitable mixing rules for application to multicomponent 
mixtures . 

The measurements focussed upon droplet combustion under 
natural and forced convection conditions. The natural convection 
tests considered high pressure droplet burning in a cold ambient 
gas environment. A wide range of alochol and n-paraffin fuels were 
tested to provide further insight into the general applicability of 
the different models of combustion and the empirical corrections for 
natural convection. 

The forced convection tests were conducted under conditions 
that were a realistic simulation of a combustion chamber environment. 
The atmospheric pressure baseline and high pressure tests provided 
burning and evaporation data for a wide range of fuels and considered 
various ambient oxygen concentrations and ambient gas temperatures. 

The tests also extended over a wide range of Reynolds numbers. The 
ambient oxygen concentration ranged from 0-13 per cent molar 
concentration and the ambient gas temperature ranged from 600°K to 
1500°K. The Reynolds number range was sufficient to provide a 
fairly stringent test of the convection correction for forced 
convection. 

The experiments indicated that methanol and ethanol were 
approaching critical combustion conditions at pressures on the order 
of 80-100 atm when burning in air under natural convection conditions. 
Both the low pressure and hi*, pressure quaternary theories predicted 
critical burning in reasonable agreement with these results. Critical 
burning conditions could not be approached for the remaining fuels 
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due to the formation of soot deposits on the sphere. It was found 
that the rate of soot formation increased with increasing ambient 
pressures and ambient gas temperatures and decreasing ambient oxygen 
concentrations for the n-paraffin fuels. 

The use of the variable property-variable specific heat, 
variable property-constant specific heat and constant property gas 
phase analyses gave essentially the same results as long as the 
respective constant properties were evaluated at average conditions 
inside and outside the flame. Parametric property variations 
caused variations in the computed results similar to those encountered 
in earlier high pressure combustion studies. [42. 55, 56] For 
porous spheres, the binary high pressure theory gave a poorer 
approximation of the quaternary high-pressure tneory, than was the 
case for high pressure droplet combustion. [42] The theoretical 
predictions for the high pressure cold gas tests were consistent 
with the trends predicted for the high pressure, high ambient gas 
temperature tests. 

Since many combustors are of the liquid spray or injection 
type, the results of this investigation will provide useful data 
concerning the effects of high pressures, temperatures, and convection 
on the phenomenon of combustion. Furthermore, the results of this 
investigation have provided fundamental burning rate data for a 
number of fuels and provided a fairly stringent test of the various 
high pressure models of combustion. 

The results of previous investigations have revealed the 
need for a more fundamental understanding of high pressure combustion 
processes. There are several interesting extensions of this 
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investigation using the present experimental set-up that can be 
considered in future studies. One of these vould be a detailed 
theoretical and experimental investigation into th * phenomenon of 
soot formation during the combustion of n-paraffin fuels at 
elevated pressures. The present experimental facilities can be 
used to make a quantitative study of the effects of ambient oxygen 
concentration, pressure, and temperature on the mechanism of fuel 
pyrolysis and soot formation. 

Another useful extension of the present work would be a 
study of the physical effects of ambient gas velocity and temperature 
on the flame stability of burning droplets at high pressures in a 
convective flow environment. Such phenomena as the transition from 
envelope to wake flames and the effect of gas velocity on the process 
in a combustion chamber environment require further investigation. 

The correlations that are presently being used for forced convection 
should be reexamined so that they may provide a better understanding 
of the physical processes that take place when a liquid fuel is 
burned in the form of droplets. 

5.2 Conclusions 

The major conclusions of this study are as follows: 

1. Both low and high pressure theories gave essentially 
the same burning rate predictions (Figures 8, 9, 10, 

22, and 23). The greatest difference between the 
theories lies in the predicted critical burning 
conditions (Tables 2 and 5). The low pressure theory 
predicts a significantly lower pressure for critical 
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combustion , at a given ambient oxygen concentration, 
than the high pressure theory. Furthermore, the low 
pressure theory predicted higher liquid surface 
temperatures at all pressures for critical combustion 
than the high pressure theory (Figures 11, 12, 24 and 
25). 

2. Discrepancies between theoretical and experimental 
burning rates were about the same as in earlier 
atmospheric pressure studies (Reference [19]). 

3. The low pressure theory gave the best temperature 
prediction for methanol, as shown in Figure 11, 
whereas the high pressure theory was best for ethanol 

and propanol-1- Over the experimental range for n-heptane 
the low pressure theory was adequate. The experimental 
pressure range for n-pentane and n-decane was not high 
enough to be conclusive (Figure 12). 

4. The experiments gave evidence of critical combustion 
for methanol and ethanol at pressures of 80-100 atm, 
in agreement with both theories. 

5. Predicted dissolved gas concentrations reached values 
as high as 60% (n-decane) for the present test 
conditions . 

6. The binary model was not a good approximation for the 
quaternary model in the present case of porous sphere 
combustion, unlike the previous findings of Lazar [43] 
for droplet combustion. 
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7. The variot i gas phase theories gave essentially the 
sane results as long as any assumed constant property 
was evaluated at an average condition. 

8. The effect of parametric property variations, as shown 
in Figure 28 and TabJe 7, was about the same as in the 
case of droplet combustion. 

9. Based on approach conditions, the multiplicative 

correction for the natural convection correction for 

the burning rate appears adequate for Grashoff numbers 
4 8 

in «:he range 10 -10 . 

10. The multiplicative correction for forced convection 
appears adequate for droplet combustion, when based 
upon approach conditions, over the Reynolds number 
range 10-800. 

11. Ambient air temperature has only a slight effect on 
the burning rates at finite ambient oxygen 
concentrations for the fuels considered in the tests 
(Figures 14 and 15) . The burning rates are very 
sensitive to temperature for very low oxygen 
concentrations in the ambient gases (Figure 19). 
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APPENDIX A 


CHECK OF ASSUMPTIONS OF THE ANALYSIS 


A. 1 Constant Total Pressure Assumption 

Reference [48] gives the following approximate expression for 
estimating the pressure changes across the boundary layer surrounding 
a spherically symmetric droplet 
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9 . 
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2RT 




DX 
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V 1 - x u 



(A. 1) 


The above equation was obtaix. d by integrating the momentum equation 
for the flow field resulting from the binary diffusion of fuel vapor 
from the droplet. In addition to the assumption of a spherically 
symmetric system, it was assumed that the molecular weight of both 
species was the same and constant fluid properties in the boundary 
layer. 

It was found in Reference [43] that for droplets as small 
as lQy in diameter at one atmosphere pressure, that Equation (A.l) 
indicates that there is less than 5% change in pressure, providing 
<_ 0.9. Since the present investigation considered much larger 
droplet sizes (0.64-1.9 cm) and pressures much greater than one 
atmosphere, the pressure variation is negligible. 

A. 2 Radiation As s umption 

A spherical model was used to estimate the contribution of 
radiant energy to the droplet from the flame. Such a model is 
obviously oversimplified since the distortion of the boundary layer 
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surrounding the droplet by convective forces will alter the assumed 
spherical symmetry of the system. Furthermore, the flames for most 
of the fuels used in the investigation were usually blue and non- 
radiant over the front half of the sphere, but yellow and highly 
radiant over the rear half. 

Following the approximation suggested by Godsave [1] the 
system consists of a hot radiating surface (the flame) surrounding 
a cooler body (the droplet). Radiation from the vapours at 
intermediate temperatures between the flame front and the drop is 
neglected in view of the sharp temperature gradient. By neglecting 
the absorption of radiation by these intermediate vapors, the two 
approximations tend to cancel each other out and the final result is 
reasonably representative of the actual process. 

For n- decane combustion at a flame temperature of 3000 C K and 
for various pressures, Reference [43] indicates that thermal 
radiation accounts for a negligible portion of the energy required to 
vaporize droplets up to lOOQy in diameter at low and moderate 
pressures, but can become appreciable at very high pressures. 

The present investigation considered the use of porous 
spheres to simulate the fuel droplets. If the radiation loss by the 
droplet is neglected, due to its relatively low temperature in 
comparison with that of the flame, the radiant energy flux absorbed 
by the droplet is: 

Q r « ^ r fc 2e f a d aT f^ (A. 2) 

In this equation e f is the emissivity of the thin zone of radiating 
gases, and a is the Stefan-Boltzmann constant. The radiating gases 
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were assumed to be carbon dioxide and water vapor. The partial 
pressures of these gases were based upon the average concentrations 
of these species at the flame zone. For porous alundum spheres, the 
absorptivity, a^, represents the combined effects of the liquid film 
and the surface absorptivities . It is known that the . Icorptivity 
of liquid films is a function of their tnickness. Owing to the 
uncertainty concerning the exact thickness of the liquid film and 
the belief that the porous sphere surface may be an important factor 
in the absorption of radiation, other investigators [7, 58, 59] 
have used an average value of 0.8 for the surface absorb tivity. 

For the case of insoluble ambient gases in the liquid phase, 
the total heat flux to the droplet for steady burning is 

Q = nL (A. 3) 

The low pressure variable property-variable specific heat 
theory was used to compute the molar burning rate and liquid surface 
temperatures. High pressure effects were considered in evaluating 
the heat of vaporization, L. For the combustion of n-heptane in 
air at a pressure of 50 atmospheres the thermal radiation was found 
to account for about 36 per cent of the energy required to vaporize 
the droplet. 
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APPENDIX B 


GAS PHASE ANALYSIS 


B. 1 Variable Property-Variable Specific Heat Analysis 

Under the assumption that compressibility effects are small 
in the gas phase, the ideal gas relationships are employed to evaluate 
the enthalpy difference terms appearing in the steady burning equations 
for conservation of energy. Substituting Equation (2.21) 

T 

B, 


(h-h £ ) = f (A ± + B T)dt - A ± (T 

i X 


- r ) + y- (T 2 -t £ 2 ) 


(B. 1) 


Substituting this expression Equation (2.8) becomes 

2 . 


n A 1 (T-T i ) + B x 
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(B. 2) 


The burning rate, n, is determined by applying Equation (2.20) 
and the boundary conditions, Equation (2.18), to Equation (B.2). Upon 
integration, the result is: 
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-2A l B l^f T^) 2 

n = (A x + B 1 T f ) (A n + B^j ’ 5 = ° (B,6) 

where 

£ = " <t> 2 = A : 2 - 2B 1 [L 1 - A x T £ - — T £ 2 ] (B. 7) 

An expression for the fuel mole fraction at the droplet 
surface may be obtained by eliminating spatial deviatives between 
Equation (2.9) and (B.2), introducing Equation (2.18), and integrating. 


The solution has three branches which are as follows: 
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(B.10) 


where 


x A - WCD) A . 


To determine the concentration of any species in terms of 
the fuel mole fraction at the droplet surface and conce. rations at 
the flame. Equation (2.9) is used. Eliminating spatial derivatives 
and introducing Equation (2.18) results in the following expression: 


hi m x if (1 " x u )# 1-2 * •••» N " 1 


(B. 11) 


In order to determine conditions at the flame surface Region 
B must be considered. The mole fraction of any species at the flame 
may be expressed in terms of the ambient oxygen concentration through 
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the use of Equation (2.17) and integrating. The result is 


a ' X if- a i 


“N 


a ' x i«r a i V a 'V 


(B. 12) 


where 


a' = f a 


(B. 13) 


i=2 

Substituting the specific heat relation Equation (2.21) 
into the energy equation leads to 
N B 

«£ °i {A i (T -V + r (l 2 - T i 2)) - + L > ■ § (B - 14) 


Letting 


a ' - L «i A i 
i=2 1 1 

(B. 15) 

N 


b' * I *i B i 

i-2 

(B. 16) 


Equation (B.14) becomes 


n[(T-Ij) «• + ^V-T^ - % + l] - 4irt 2 X B g 


(B. 17) 


An expression for the combustion temperature in terms of the 
ambient oxygen concentration may be obtained by eliminating spatial 
derivatives between the conservation of energy and conservation of 
species equation for Region B and integrating between the limits 


T - T f 
T - T 

00 



This results in the following expressions: 


(B. 18) 



I 
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-■ft—,/ 

I 


/ \ 
(b'T. + a' + Y)(b'T f + a’ - y) 


X B / Yet' 


(b'T. + a' - y)(b'T + a' + y) 




y > 0 


(B. 19) 


2a' X 




B I -1 

exp < tan 


a' + b'T 


f J . -! * a ' +b ’ T - 




- tan 


r 


In 


% 


2a'x B b , (T f -T o ) 


V^V/ <•' + b'T f )(a' + b'T.) ’ 


where 


in 


% 


°y 


-a'XN 


< t >’ 2 » 0 


(j)' 2 > 0 (B.20) 

(B. 21) 


Y - " <J>' 2 * a' 2 - 2b' (1^ - Q z - a’T z - y^ 2 ) 


(B. 22) 


An expression for the combustion radius can be obtained 
through the use of Equations (B.17) and (2.19) and integrating. 

There are again three branches to the solution depending on 
the value of y. They assume the form 


4irX 


BH 


r f 


in 


< T «-V a ' + b'/2(T ao 2 -T Jl 2 ) - Q l + ^ 

( W a ' + r (I f 2 - T t 2) - «i + h 


(3.23) 


, f(a' +b'T f - y)(a' + b'T. + y)] 


n ’ “ y ln [ (a* + b'T f + y)(a' + b'T. - y) J 


n' 


2a' r -1 
jr- I tan 


(a' + b'T.) , (a' + b'T ) 

f . -1 00 

tar 


l ’ ■ 


r 

2a'b’(T f -T .) 


r 


, y > 0 (B. 24) 

, <fr’ 2 > 0 (B.25) 


(a' + b'T f )(a' + b'T.) » 


2 

y - 0 


(B.26) 
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B.2 Variable Property-Constant Specific Heat Analysis 

If is set equal to zero in Equation (2.21) the solution 
corresponds to the constant specific heat version of the gas phase 
analysis. For this case the energy equation for Region A becomes 


n[A 1 (T-T Jl ) + L 2 ] 


X 

7j^T4Trr 

l Z 


2 dT 
dr 


(B. 27) 


Integrating this expression between and r for r < r^ leads to 
the following result for the combustion radius: 


°Vi 


r £ 
1 - — 

"f 


WV ~ (L r A iV in {1 + l (T f~ T £ )} (B - 28) 


An expression for the liquid surface mole fraction may be 
obtained by eliminating spatial derivatives between the equations of 
conservation of energy and species and integrating between the limits 
given in Equation (2.18). The result is 




1 - 


-> 1/A 1 X A 


h + VW 


(B. 29) 


Equation (B.ll) is unchanged for this case and remains 


X u - X lf (l-X u ) 1*2, .... N-l 


Considering the outer region. Equation (B.17) with b'*0 

becomes 


n[a' (T-T^) - (fy + I^J 


47rr 2 X 


T_ dT 
BZ T^ dr 


(B.30) 


When this equation is Integrated between the limits 
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When this equation is integrated between the limits 


r " r f T " T f 


r * r 


T - T 


(B. 31) 


the resulting expression for the burning rate is 


4l,r £ V 


EW - <h - «i 


a’Tj) In 


[»’<W - <U + L ll 

a' (I f -T t ) - Qj + 


(B. 32) 


An expression for the combustion temperature may be developed 
through the use of the conservation of energy and species equations. 
The result is 



% 

a N _X N“ C ‘ . 



.»'<w - «i- v 


(B. 33) 


B. 3 Constant Property Analysis 

In the assumptions that the fuel thermal conductivity and 
specific heat are independent of temperature are made along with the 
assumptions given in Chapter II, the analysis will correspond 
essentially to that of Reference [1], With the exception of the 
combustion temperature, expressions for the unknowns must be rederived 
from the steady burning equations presented in Chapter II. Constant 
average specific heats are used for the fuel vapor, oxidizer and 
products, designated as Cpp, Cp Q , and Cp^, respectively. The 
derivation of the equations is developed in the same manner as 
presented in Section B.l and the details are omitted nere. 
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The resulting expressions for liquid phase mole fractions, 
combustion temperature, flame radius, and burning rate are: 




1-<1 + 


Cp F ( W 


-A/CDCp, 




X 2* = X “ X U 


(B. 34) 


(B. 35) 


in 


»’<W -«* + L i 

al 


“N J 

< W - + \ 

a' 

* 

1 

p 

L_l 


47iAb 


•'<VV - Q * + L i 


r f Cp Q n 
— - — 1-rrf- - pi-in fl 

r SL r f J 4 X A Cp F l 


W- Qi + "V 

^VVVj 


(B. 36) 


(B. 37) 


(B. 38) 


APPENDIX C 


PHASE EQUILIBRIUM RELATIONS 


C . I Mixing Rules 

The Redlich-Kwong equation of state is 


^ RT 
P - TT-To 


V-b x°' 5 V(V+b 9 ) 


(C.l) 


wuicn is cubic in V. The largest positive root is used in 
calculations . 

The mixing rules of Prausnitz and Chueh [51] are used to 
apply the above equation to mixtures. For a mixture of N components: 

N N 


where 


and 


- y y x.x.a* 




“ai R \l 

P ci 


2.5 




2 2 5 

(fi , + u ,)r t ., J 

a* aj cij 
cij 


(C. 2) 


(C. 3) 


(C. 4) 


X i b i 


. W 


ci 


Z ,. RT .. 

P - .-g.U _ c il 

ci J V cij 


(C.5) 


(C.6) 


(C. 7) 


V ctJ • <\l + V CJ )/2 


(C.8) 




, . = 0.291 - 0.04(o)j + to . ) 
-ij 1 J 


(C.9) 


T , - /T .T .'(l-k. ,) 

clj cl CJ ij 


(C.10) 


C . 2 Component Fugacities 

Component fugacities were determined by substituting 
Equation (C.l) and the mixing rules, Equations (C.2) to (C.10), into 
Equation (2.54) and integrating. The result for a component i in a 
mixture of N components is given by the following expression: 


I v \ b l 2 X ^ a j i 

£n(f ± /X ± P) - ^nj v _ b 6 | + — b - - RT 3/2 b . 


5 V+b«\l, a ° b I 

I U RT 3 / 2 b « 2 


o PV 
lu RT 


(C.ll) 


The molar volume, N, is determined from the equation of state which 
is cubic in V. 


heat of Vaporization 


The heat of /aporization of a component i in a mixture, 
is the difference between the partial molar enthalpy in the vapor and 


liquid phases, i.e.. 


l i • s i v - 


(C. 12) 


The partial solar anthalpy of component i in a mixture is determined 
by the folloving thermodynamic relation: 


!Lii. i_L II. 
ri 2 21 r v 


(C. 13) 


In the above expression h* represents the partial molar enthalpy in 
the ideal state. When Fquation (C.12) is differentiated and the 
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resulting expression substituted into Equation (C .13) the following 
expression results: 


h » “ h C A - 

= (A-l) \ + Y + B 

RT 


(C.14) 


where 


N 

2 V a! # X. o K o 

_ -b° _ b i V + 1=1 l _ j_ 

v " b ° (V-b ®) 2 RT 1 /, 5 (V+b°) RT 1 * 5 (V+b 5 )" 


(C. 15) 


B = 


2RT 2 'V 


2 T X.a! 

,4ri J J 1 
. lJ" 1 


0 i a b 

ii £n l V b 


_ Oi o 

i f I V+b 0 i 

Zn '~v~ + V+bT 


(3Z 2 - 2Z + D)C = 2jr- + 

(RT) ‘ 


2T 


3/2 


+ b°R 


3 . 5a®b°P 2 

rV - 5 


(C. 16) 


(C. 17) 


D = 


.1/2 


- b’RT - b Pi 


(RT)' 


(C. 18) 


and 


PV 

7 = — — 

4 RT 


(C. 19) 


is the mixture compressibility tactor. 

By solving Equation (C.14) for each phase and substituting 
the results into Equation (C.12), the heat of vaporization for each 
component in the mixture is determined* 


C.4 Physical Constants 

The pure component constants P c> T c> V^, and to employed in 

the calculations are given in Table 8. 

The binary interaction parameters, k^, for the alcohol 
and n-paraffin fuels used in the investigation are presented in Table 9 




I 

i 
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Table 8 


Pure Component Constants 


Component 


P (ATM) 
c 


CC 


gm-molel 


T (°K) 
c 


00 


M 


. . E EL— 

gm-mole 


N 2 

33.54 

C °2 

72.80 

h 2 o 

218.3 

CH 3 OH 

78.5 

C 2 H 5 OH 

63.0 

C 3 H ? 0H 

51.0 

C 5 H 12 

33.3 

C 7 H 16 

27.0 

C 10 H 22 

20.8 


90.1 

126.26 

94.0 

304.16 

55.2 

647.4 

118.0 

513.2 

161.3 

516.0 

220.0 

540.7 

311.0 

469.5 

431.9 

540.2 

620.0 

619.0 


.040 

28.016 

.225 

44.01 

.344 

18.02 

.556 

32.04 

.635 

46.07 

.600 

60.09 

.252 

72.15 

.349 

100.20 

.479 

142.28 










Binary Interaction Parameters 


Component 

N 2 

co 2 

H 2 0 

Methanol 

.10 

.08 

.15 

Ethanol 

.15 

.11 

.20 

Propanol-1 

.20 

.16 

.25 

N-Pentane 

.25 

.20 

.30 

N-He tane 

.35 

.29 

.40 

N-Decane 

.50 

.43 

.55 

Note: k.. = k.. and k . , = 0 

ij ji ii 


i 





APPENDIX D 


PHYSICAL PROPERTIES 

D. 1 References for Physical Properties 

The sources of property data and correlations were the same 
as in earlier studies. [42, 43, 56] These sources are listed in 
Table 10. The specific correlations are given in the following 
sections. 

For the low pressure theory calculations the combustion 
products were treated as a single species. The properties for this 
effective species was obtained by averaging the properties for an 
equimolar mixture of carbon dioxide, nitrogen and water vapor. 

In calculations that required constant average values, the 
average temperature was used in the particular region. For Regions 
A and B these correlations assume the form 

T A = (T £ + T f )/2 (D,1) 

and 

T fi - (T f + TJ/2 (D, 2) 

D. 2 Liquid Phase Properties 

The liquid molar density, vapor pressure, and heat of 
vaporization employed in the low pressure calculations were correlated 
with the following equations: 

C J1 • (C 1 - C 2 T H )/M l* 5 ^) (D - 3) 

The vapor pressure correlation for the alcohols assumed the form 

P v - exp (C j - Cj/Tj) 


D.4) 
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Table 10 


References for the Physical Properties 




References 


Property 

Paraffins 

Alochols 

Ambient Gas 


61 

62 

— 

C 

61 

52 a 

63 

p 

D 

52 b 

52 b 

52 b 

L 

52 C 

52 C 

— 

P 

61 d 

64 

mm — 

V 

X 

60 

63 

63 

u 

60 

63 

63 


Computed, Rihani and Doraiswamy Method 
^Computed Fuller, Schettler, Giddings Method 

Q 

Watson Correlation 
^Antoine Correlation 
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and for the n-paraffins 


P v = exp{C 5 - C 6 /(T a - C 7 )} (D. 5) 

The heat of vaporization correlation was 

L = C g (l - V T c )0 * 38 (D.6) 

The constants to Cg are tabulated in Table 11. The 
molecular weight, M, and critical temperature, T , for each of the 
components considered in the calculations are also listed. In the 
above calculations the liquid temperature, T^, is in degrees Kelvin. 

D. 3 Gas Phase Properties 

The specific heat of each species in the gas phase was 
represented by the equation: 

C p = A + BT (cal/m-mole °K) (D.7) 

The specific heat constants, A and B, for each species are 
tabulated in Table 12. 

The thermal conductivity correlations employed for Regions 
A and B are as follows: 


X A ■ \U <T/ V (D ' 8) 

and 

X B - X BJ2 (T/ V (D,9) 

The constants X ^ and X g ^ represent the thermal conductivities of the 

gas mixture, evaluated at the temperature, T ^ , for Regions A and B, 

avg 

respectively. The correlation used assumed the 



Constants in the Liquid Phase Property Equations 
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Table 12 


Specific Heat and Thermal Conductivity Constants 


Component 

A 

B 

C 9 

C 10 x 10 

C^ x 10 

N-Pentane 

31.75 

28.86 

— 

— 

215.00 

N-Heptane 

44.10 

40.00 

— 

— 

205.00 

N-Decane 

62.70 

56.90 

— 

— 

190.00 

Methanol 

11.00 

9.00 

— 

— 

340.00 

Ethanol 

17.60 

13.20 

— 

— 

305.00 

Propanol-1 

22.00 

19.20 

— 

— 

305.00 

Water Vapor 

7.90 

1.80 

-47 . 40 

332.50 

— 

Carbon Dioxide 

9.00 

2.50 

11.90 

139.70 

— 

Nitrogen 

6.90 

.75 

31.80 

123.20 

— 

Oxygen 

7.10 

1.00 

45.50 

132.20 

— 


, -fr 
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N X. X . 

Z i i 

. N 


(D. 10) 


T x. e, . 

A J lj 


where 


..i44f /2 U*f&r 

1 JS\ V l 'V Iv 


(D.ll) 


and 


X^ = mole fraction of i 
= molecular weight of i 
y^ ■ viscosity of i 

The species thermal conductivities, X^, were computed from the 


correlation 


X x 10 6 - C 9 + C 1£J x 10“ 3 T 


(D. 12) 


for the non-fuel species and from the correlation 

X X 10 6 * C U (T - 140) 


(D. 13) 


foi the fuel species. 

The constants to are tabulated in Table 12 for all of 
the fuels considered in the investigation as well as the values for 
oxygen and the combustion products. 

The gas phase molar density was computed from the ideal gas 
equation of state 


C » P/RT (D.14) 

The values for the viscosity of each species is given by 


the correlation 


m 


y i = 10 6(C 12 + C 13 x 10 3t) gm/cm-sec (D.15) 


The constants used in the above expression are listed in 

Table 13. 

The binary diffusion coefficients for Regions A and B were 
computed according to the Fuller, Schettler, and Giddings [52] 
empirical correlation. This expression assumes the form 


10" 3 T 1,75 (1/M i + 1/Mj) 1/2 
ij P[(^v) 1 1/3 + (Ev)j 1/3 ] 2 


(D. 16) 


with D in cm /sec and P and T in atm and °K, respectively. The 
values of the sums of atomic diffusion volumes for each of the gas 
species are listed in Table 13. The mixture dif fusivities were 
obtained from the correlation 


^i(mean) 


(1-X.) 


X, 

ij 


w . 


(D.17) 


D.4 


Heat of Reaction 


The standard heat of reaction for the six fuels used in 
the tests are tabulated in Table 14. The reference temperature was 
taken as 298°K. 

D.5 Ambient Gas Properties 

The Reynolds number and Prandtl number appearing in the 
convection correction, Equations (2.50) and (2.51) were evaluated 
for the ambient gas mixture. The properties appearing in these 
two dimensionless numbers were determined for a gas mixture of N 
component as follows: 
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Table 13 


Viscosity Constants and Atomic Diffusion Volumes 



Component 

C 12 

C 13 

(Zvp 


Methanol 

65.00 

235.00 

29.90 

Ethanol 

60.00 

190.00 

50.36 

Propanol 

60.00 

190.00 

70.80 

Pentane 

80.00 

125.00 

106.26 

Heptane 

55.00 

120.00 

147.18 

Decane 

35.00 

50.00 

208.56 

Carbon Dioxide 

152.00 

235.00 

26.90 

Water Vapor 

136.00 

260.00 

12.70 

Nitrogen 

152.00 

235.00 

17.90 

Oxygen 

202.00 

273,00 

16.60 
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Table 14 


Heats of Reaction 




Component 

Standard Heat of Reaction 
(kcal/g-mole) 


N-Pentane 

782.0 

N- Heptane 

1075.9 

N-Decane 

1516.6 

Methanol 

161.7 

Ethanol 

305.5 

Propanol-1 

452.2 



Specific Heat 


L Vi c pi'' X, Vi 


(D. 20) 


C pi " °19 + C 20 T » (cal/ 8 rc ° K) 


(D. 21) 


Viscosity 


w ' L V 


(D. 22) 


w i “ C 21 + C 22 T » 


(D. 23) 


Thermal Conductivity 


X * l 


(D. 24) 


X i ' C 23 + C 2A T - 


(D. 25) 


Density 


P - P t X.M JU (gm/cc) 

i-1 1 1 


(D.26) 


The constants C,_ to C-, used in the above correlations are 
19 24 


tabulated in Table IS. 


s»-- * m 




APPENDIX E 


EXPERIMENTAL RESULTS 


Table 16 


Experimental Burning Rates for Natural Convection Tests with Fuel 
Inlet and Ambient Air Temperature of 300 °K 


Fuel 

D S 

(cm) 

P 

(atm) 

mp x 10^ 
(gm/sec) 


Methanol 

0. 

.64 

1.00 

0.72 

A: 

Lr 




2.18 

0.94 






10.20 

1.43 






17.35 

1.99 




1 


23.81 

2.09 






28.23 

2.14 


t 

i 




36.74 

2.33 


i. 


> 

f 

37.42 

2.38 

4 , 

Methanol 

0, 

.95 

1.00 

0.82 

Air 




4.08 

2.38 






8.84 

3.10 






16.24 

3.45 






29.39 

4.26 






32.65 

4.91 






42.18 

4.76 






44.22 

5.67 






61.22 

5.74 




> 

if 

78.23 

7.19 

1 

Methanol 

1. 

,90 

1.00 

2.93 

Air 


1.70 

2.39 

2.72 

4.08 

6.80 

12.25 

14.97 

21.09 

27.21 

28.23 

34.01 


3.99 

4.51 

4.49 

5.41 

6.20 

7.92 

8.62 

9.90 

10.56 

11.48 

12.74 
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Table 16 (Continued) 



Fuel 

] 


P 

2 

m F x 10 

X o 


(cm) 

(atm) 

(gm/sec) 

°2°° 

Ethanol 

o 

.95 

l.Ou 

1.10 

Air 




1.76 

1.48 






3.40 

1.88 






4.08 

2.28 






4.76 

2.73 






6.80 

3.15 






8.11 

3.28 






16.32 

3.62 






33.33 

7.17 






41.49 

8.31 




\ 

/ 

71.42 

9.86 

\ 

/ 

Propanol-1 

0. 

.95 

1.00 

1.34 

Air 




1.36 

1 42 






2.17 

1.69 






3.40 

1.96 






4.08 

2.03 






5.44 

2.11 






8.02 

2.38 






11.56 

2.57 






12.92 

2.59 






20.40 

3.57 






27.21 

3.82 




\ 

/ 

31.97 

4.72 

\ 

/ 

N-Pentane 

0, 

.95 

1.00 

1.83 

Air 



! 

1.29 

2.22 






2. 72 

3.27 






6.80 

3.68 






13.60 

4.42 






27.21 

4.98 




\ 

/ 

51.02 

5.64 

\ 

t 

N-Heptane 

0, 

► 64 

1.00 

0.73 

Air 




1.50 

0.74 






2.85 

1.06 






5.60 

1.35 






8.20 

1.45 






10.50 

1.49 






14.39 

2.52 




N 

/ 

21.77 

2.94 

\ 

/ 



138 


Table 16 (Continued) 


N-Heptane 


N- heptane 


N-Decane 


N-Decane 


N-Decane 


p 

(atm) 

mp x 10^ 
(gm/sec) 

1.09 

1.56 

2.04 

1.69 

4.08 

2.06 

8.16 

2.66 

10.88 

2.71 

18.71 

3.04 

27.21 

3.57 

04.42 

3.79 

1.00 

3.94 

1.91 

4.53 

2.38 

4.59 

3.06 

5.43 

4.08 

5.66 

6.12 

6.46 

8.16 

8.09 

9.52 

8.29 

12.93 

10.65 

17.00 

13.16 

1.00 

0.49 

1.35 

0.57 

4.08 

0.98 

5 . t '<■ 

1.35 

8.33 

1.38 

18.37 

1.47 

1.09 

1.21 

1.36 

1.41 

2.72 

1.96 

6.26 

2.17 

13.61 

3.29 

32.65 

3.98 

1.00 

3.23 

2.04 

4.24 

3.40 

4.64 

5.98 

5.38 

7.48 

6.05 

8.03 

6.32 

8.84 

6.32 

9.52 

6.92 

12.24 

7.26 

12.93 

7.64 
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Table 17 


Experimental Liquid Surface Temperatures for Natural Convection 
Tests with Fuel Inlet and Ambient Air Temperature of 300°K 


Fuel 


D 

(atm) 


T£ 

(°K) 


2 “ 




Methanol 

1.00 

328 

Air 




2.04 

343 






6.80 

379 






8.16 

387 






10.88 

397 






12.24 

402 






14.96 

411 






21.76 

432 



; 



38.09 

446 



*■ 

t 


44.89 

460 



( 



47.61 

465 



i . 



61.22 

486 



t 

4 


78.23 

497 

\ 

/ 

t I 

i 

Ethanol 

1.00 

341 

A: 

.r 

i 

\ 


2.72 

359 



4 

v i 


4.08 

372 



r f 

i* 


5.44 

384 



i 



6.80 

391 



* 

f 

1 

a 

f 

4 


8.16 

397 



r 

t 


10.84 

407 



i 

i 


13.60 

410 



i 

i 


18.70 

415 




i 


34.01 

426 



c 

? 



48.97 

442 



* 

f 



61.22 

453 

\ 

r 

- 


Propanol-1 

1.49 

342 

Aj 

.r 

i 



2.04 

367 



« , 



3.53 

370 






3.80 

389 



* • t ■ 



6.80 

397 



* 



7.48 

409 

i 





20.40 

447 






29.53 

461 



t 



43.53 

468 

\ 

t 













^ -! '"I? . 






r 






S’ 

f. ■ 
t 

s * 

f 

i 

* 

i 















1 

Fuel 

P 

T£ 

v ; 


(atm) 

( C K) 

N-Pentane 

1.00 

298 

1 

Air 1 


2.38 

337 


1 


3.40 

339 




6.12 

341 




6.80 

342 




9.52 

353 




11.22 

356 




13.60 

wz 


1 


20.40 

413 

N 

t 

N-Heptane 

1.00 

360 

Air 

2.38 

378 




3.06 

383 




5.10 

422 




6.12 

436 




7.48 

444 




8.84 

433 


\ 


11.56 

448 


t 

| 


12.92 

453 




13.60 

454 




15.64 

475 




17.00 

469 

\ 

! 

N-Decane 

1.09 

432 

Air 


1.49 

435 




2.38 

453 




3.40 471 

6 . 80 480 

7.48 493 

9.52 514 

15.75 540 

18.72 561 

24.28 565 
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Table 18 


Experimental Vaporization Rates for Forced Convection 
Tests at Atmospheric Pressure 


Fuel 


T 

(°K) 

Reynolds 

Number 

in 2 

mp x 10 
(gm/sec) 

v 

610 

116 

0.32 

Evap . 

630 

131 

0.36 

I 


790 

101 

0.46 



795 

142 

0.31 



800 

101 

0.43 



845 

174 

0.66 



895 

152 

0.73 



1090 

52 

0.85 



1090 

116 

1.26 



1130 

112 

1.38 



1160 

49 

0.99 

1 

i 


1180 

64 

1.20 

i 

i 


1255 

119 

1.87 



1270 

46 

1.46 



1320 

75 

1.54 



1350 

141 

2.40 



1480 

171 

2.71 

\ 

f 

1150 

110 

2.13 

0.095 

1160 

55 

1.83 

1 


1310 

128 

2.29 



1340 

100 

2.17 


t 

780 

106 

1.10 

Evap . 

795 

75 

1.08 



1000 

114 

1.00 



1060 

110 

1.12 



1230 

99 

1.63 



1365 

112 

2.06 

> 


1100 

128 

1.26 

0.095 

1130 

166 

1.66 



1265 

78 

2.04 



1270 

116 

1.97 

> 

/ 

610 

119 

0.78 

Evap. 

615 

155 

0.99 



1045 

52 

0.62 



1090 

55 

1.26 



1145 

82 

0.83 

| 

\ 

/ 


Methanol 


Methanol 


Ethanol 


Ethanol 


Propanol- 1 
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Table 18 (Continued) 



Fuel 

T 

_ oo 

Reynolds 

2 

Ap X 10 

X n 



( K) 

Number 

(gm/sec) 

2“ 

Propanol-1 

1200 

79 

0 . 8i. 

Evap . 


1265 

50 

1.52 




1285 

57 

1.04 




1310 

115 

2.19 




1365 

87 

1.60 




1395 

92 

2.10 




1530 

75 

2.87 

\ 

/ 

Propanol-1 

1215 

67 

1.71 

0.095 


1255 

54 

2.17 




1365 

74 

2.50 




1380 

97 

2.62 




1380 

93 

2.70 




1450 

182 

3.55 

N 


N-Pentane 

980 

118 

1.00 

Evap . 


1090 

109 

1.95 




1255 

99 

2.31 




1035 

134 

2.23 




1090 

115 

2.95 




1160 

111 

2.63 




1255 

118 

3.55 

\ 

/ 

N-Heptane 

885 

71 

0.92 

Evap . 


910 

137 

1.16 




1155 

87 

1.06 




1200 

99 

1.46 




1255 

90 

1.81 




1255 

115 

2.10 

\ 

f 

N-Heptane 

1090 

111 

1.69 

0.095 


1170 

134 

2.72 




1245 

119 

2.86 




1255 

92 

2.64 




1255 

79 

1.93 

N 

/ 

N-Decane 

910 

75 

0.61 

Evap. 


920 

76 

0.54 




1000 

152 

1.19 




1030 

112 

1.29 




1255 

79 

1.03 




1330 

104 

1.49 




1365 

143 

2.03 

> 

/ 



Table 18 (Continued) 
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* 

v, 


t 


t 

i 



Fuel 

T 

Reynolds 

mp x 10^ 

X 0 


(°K) 

Number 

(gm/sec) 

°2“ 

N-Decane 

1130 

126 

1.42 

0.095 


1220 

80 

1.62 



1310 

115 

1.59 

N/ 
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Table 19 


Experimental Liquid Surface Temperatures for Forced Convection 
Tests at Atmospheric Pressure 


Fuel 

T 

CO 

(°K) 

Reynolds 

Number 

To 

(°K> 

Ethanol 

1000 

114 

335 


1060 

110 

335 


1230 

99 

341 


1365 

112 

342 

Ethanol 

1100 

128 

343 


1130 

166 

343 


1265 

78 

344 


1270 

116 

345 

N-Penta.ie 

795 

160 

293 


820 

160 

293 


920 

107 

294 


1090 

130 

299 


1255 

117 

303 


1310 

117 

306 

N-Pentane 

1085 

86 

305 


1130 

127 

306 


1160 

111 

306 


2 “ 

Evap. 


0.095 


0.095 
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Table 20 


High Pressure Experimental Vaporization Rates and Liquid 
Surface Temperatures for Forced Convection Tests 








Fuel 

P 

(atm) 

Reynolds 

Number 

TO 2 

mp x 10 
(gm/sec) 

T£ 

(°K) 

3 

CM 

O 

X 


T - 1255°K 
00 


Ethanol 


I 


Ethanol 


t 

i 

f 


Ethanol 


1.00 

99 

1.63 

335 

Evap . 

2.99 



361 



6.80 

167 

2.16 

381 



6.80 

167 

2.07 

381 



13.61 

312 

2.36 

395 



20.40 

315 

2.74 

424 



27.20 

406 

2.54 

443 



30.61 

416 


440 

\ 

/ 


T =■ 
00 

1145°K 




1.00 

165 

1.66 

343 

0.095 

1.36 

82 

1.75 

347 



4.76 

145 

1.53 

381 



4.76 

145 

1.67 

361 



6.80 

243 

2.39 

391 



13.61 

325 

3.03 

399 



20.40 

461 

3.15 

417 



27.20 

617 

3.92 

439 

\ 


1.00 

156 

2.05 

336 

o.: 

.30 

7.10 

276 

3.00 

390 



13.00 

437 

3.43 

418 



20.40 

657 

4.50 

430 

> 

/ 

1.00 

129 

1.06 

350 

Evap. 

1.63 

188 

2.12 




2.72 

66 

1.86 

388 



2.72 

188 

2.03 




6.80 

236 


404 



6.80 

236 

2.00 




13.70 

236 

2.30 




13.70 

314 

2.52 




14.28 

314 

2.30 

441 



20.40 

448 

2.61 




23.20 

534 

3.24 




27.55 

640 

2.36 




34.01 

672 

3.24 

486 

> 

/ 


N-Heptane 
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Table 

20 (Continued) 





Fuel 

p 

Reynolds 

ntp x 10^ 

T* 

x o 



(atm) 

Number 

(gm/sec) 

(°K) 

l 

>00 



T = 

1145°K 






00 





N-Keptane 

1.00 

101 

1.77 

364 

o.c 

)95 


2.12 

188 

2.12 

372 




3.40 

157 

3.02 

390 




5.44 

251 

2.72 

413 


l 


6.80 

236 

2.03 





13.61 

314 

3.16 

457 


1 


14.28 

314 

3.16 

459 




T = 

1145°K 






00 





N-Heptane 

1.00 

150 

2.52 


0.] 

L30 


1.00 

134 

1.53 





2.04 

178 

2.07 

373 




7.48 

267 

2.19 

413 




8.75 

276 

2.26 

413 

N 

/ 


i 

i 

v 



I 



